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1 . 0 GRANT OBJECTIVES 


The primary objective of this grant has been to compare the general 
circulation of the Northern and Southern Hemispheres using Nimbus satellite 
data. During the past year Dr. Robert F. Adler has completed his Ph.D. 
degree working on this problem. His dissertation, entitled "A Comparison of 
the Structure and Flow Characteristics of the Upper Troposphere and Strato- 
sphere of the Northern and Southern Hemispheres” is enclosed in the Appendix 
of this report. 
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2.0 DISCUSSION OF RESULTS AND CONCLUSIONS 

2.1 Interhemi spheric Comparison of Atmospheric Circulation 
Features as Evaluated from Nimbus 3 SIRS Data. 

The results summarized here are presented in detail in the 

Appendix. 

The basic purpose of Dr. Adler’s study was to compare the structure and 
flow characteristics of the upper troposphere and stratosphere of the 
Northern and S'outhern Hemispheres. Atmospheric structure obtained from 
satellite, multi-channel radiance data has provided the primary basis 
for the comparisons. The results of this study confirm some previously- 
made conclusions and assert new conclusions regarding inter-hemispheric 
comparisons and the Southern Hemisphere circulation. 

The satellite, multi-channel radiance data have been shown to be 
a valuable asset in the study of the general circulation. In this 
study they have provided information on atmospheric structure with 
equal, or nearly equal, time and space resolutions in the two hemi- 
spheres. This investigation is, therefore, the first comparison of 
circulations of the Northern and Southern Hemispheres based on a data 
system having comparable resolution and quality of data above and below 
the equator. The radiosonde network of the Southern Hemisphere has 
nearly an order of magnitude fewer stations than the Northern Hemisphere 
version. Therefore, the comparison between the hemispheres of statis- 
tics based on these radiosonde systems may be open to question. This 
may be especially true of quantities whose estimation requires good 
horizontal resolution on the time scale of a day. 
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Most previous observational studies of the atmosphere's general 
circulation have centered on the Northern Hemisphere, because of the 
better upper-air network there. However, in order to understand the 
general circulation on a global basis we must first develop a knowledge 
of the characteristics of the circulation in both hemispheres. Because 
there is little hope that the number of Southern Hemisphere radiosonde 
stations will increase to become comparable to the network in the 
Northern Hemisphere, there will be continued heavy reliance on satellite 
data to study the general circulation below the equator. Numerical 
models of the global general circulation in the troposphere and strato- 
sphere now exist. In the Northern Hemisphere these models can be 
verified against statistics compiled from the relatively dense radio- 
sonde network; in the Southern Hemisphere this verification will have to 
depend more heavily on satellite-based studies, especially for eddy 
terms, and especially for layers above the mid-troposphere. Hopefully, 
the current investigation has provided a beginning in that direction. 

The results of this investigation do show that the satellite 
radiance data, as used here, have two main limitations in this type of 
study. First, because of cloud effects in the lower levels, results 
for the most part are restricted to layers above 500 mb. Even above 
500 mb results improve with height into the stratosphere. The second 
limitation is the tendency of the satellite-derived atmospheric structure 
to underestimate the amplitude of eddies. 

The major results and conclusions of this study are summarized in 
the following paragraphs. 

Zonal ly- averaged vertical motion fields in the stratosphere calcu- 
lated in both hemispheres confirm a two-cell structure to the mean 
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meridional circulation in Northern Hemisphere winter with rising motion 

over the North Pole. In the Southern Hemisphere winter, the calcula- 

) 

tions also show a two-cell structure in the mid-stratosphere, but with 
the main ascending branch located at approximately 70°S. The mean 
circulation in the vicinity of the Antarctic stratospheric jet stream 
is indirect, thus the maintenance of the jet stream must depend on the 
large-scale eddies. The calculation of the mean meridional circulation 
of the Southern Hemisphere stratosphere presented in this paper is the 
first realistic estimation of that circulation. That is, it is the 
first to include the very important influence of large-scale eddies. 

There are significant variations during the Southern Hemisphere 
winter in relation to stratospheric synoptic events. During a minor 
midwinter warming, the flow is dominated by a stationary wave number 1 
pattern, the eddy flux of heat at 60°S is large, and the vertical 
motion magnitudes are also large. During a "quiet" period, the flow is 
dominated by a transient wave number 2 pattern, the eddy flux is 
relatively small, and a descending branch of the mean meridional circu- 
lation appears over the South Pole. Consequently, during this quiet 
period there appears to be three instead of two cells in the mean 
meridional circulation. 

A previous calculation based on IGY data in the Southern Hemisphere 
concluded that the Southern Hemisphere mid- latitude eddy flux of heat 
is approximately one half as large as that in the Northern Hemisphere. 
Another study, based on energy budget considerations, concludes that 
the two values are about the same. The results of the present study 
indicate that the eddy flux is much closer to being equal in the two 
hemispheres than to being different by a factor of two. 



5 


The existence of longitudinal variations in the pattern of eddy 
heat flux in southern mid- latitudes, as there are in northern mid- 
latitudes, is also shown. These variations are best defined in both 
hemispheres during winter. The Northern Hemisphere variations are 
closely tied to standing waves; in the Southern Hemisphere the varia- 
tions are due more to longitudinal variations of the strength of 
transient eddies. In winter, at 45°S, three peaks are well defined, 
seemingly associated with the three land areas, Australia, South 
America and Africa. 

The relative importance of standing eddies in the two hemispheres 
has been examined closely in this study. The calculations confirm that, 
on a hemispheric-averaged basis, the Southern Hemisphere has a lower 
percentage of AE and KE contained in standing eddies. In the Southern 
Hemisphere the percentage of energy contained in the standing eddies is 
smaller in summer than in winter; in the Northern Hemisphere, the 
percentage either remains the same or increases from winter to summer, 
depending on the layer and parameter under consideration. This differ- 
ence reflects the dominance of the Asian summer monsoon in the Northern 
Hemisphere. 

The distribution with latitude in the upper troposphere of the 
ratio of standing eddy temperature variance to the total variance shows 
large inter-hemispheric differences. This parameter indicates the 
importance of standing eddies as a function of latitude. In winter, 
poleward of 30°, the Northern Hemisphere shows one peak at about 50°N; 
the Southern Hemisphere has two peaks, one at 40°S and one over 
Antarctica. In mid-latitudes (40-60°), the Southern Hemisphere value is 
lower than that in the Northern Hemisphere; in high latitudes (>60°), 
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the reverse is true. In summer, in low latitudes (<35°) , both hemi- 
spheres show maxima associated with the summer monsoon in each 
hemisphere. The Northern Hemisphere maximum is much larger than the 
Southern Hemisphere one. At all latitudes in summer the ratio is 
larger in the Northern Hemisphere. 

In winter, at the latitude of maximum eddy kinetic energy (~30°), 
the standing eddies contribute substantially to the total. In the 
Southern Hemisphere, the percentage contribution is about 30-6; in the 
Northern Hemisphere it is about 40%. These percentages are larger than 
the hemispheric-averaged values. In summer the percentages are less 
in both hemispheres . 

In the Northern Hemisphere during winter, about 50% of the eddy 
heat transport in mid-latitudes is accomplished by the standing eddies. 
In the Southern Hemisphere, the percentage contribution is much less. 
According to calculations in this study, the Southern Hemisphere per- 
centage is about 20%. Other previous studies have indicated a lower 
percentage . 

.The satellite-based zonal mean temperature and mass structure 
con&jrm the inter-hemispheric differences previously noted from 
climatological information. The upper tropospheric temperature gradient 
between 20° and 80° latitude is larger in the Southern Hemisphere in 
both winter and summer, with the largest disparity occurring in summer. 
The same difference is evident in the distribution of 200 mb geopoten- 
tial height. These variations are also reflected in the computed values 
of zonal available potential energy (AZ) and zonal kinetic energy (KZ) . 
In summer the AZ of the Southern Hemisphere is twice as large as that 
* in the Northern Hemisphere. Because the winter difference is relatively 
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small, there is an annual cycle of global AZ with peak values occurring 
during the northern winter. 

The contributions from the upper tropospheric and tropopause layers 
to the eddy available potential energy (AE) and eddy kinetic energy (KE) 
do not show the large inter-hemispheric differences that are evident in 
the -zonal components, AZ and KZ. Therefore, the ratios of AZ to AE and 
KZ to KE are larger in the Southern Hemisphere. 

The winter stratospheres of the Northern and Southern Hemispheres 
are also studied in relation to the difference in breakdown clima- 
tologies. The temperature structure of the winter stratospheric vortices 
are shown to vary between the hemispheres in a way that a previously- 
deriyed instability criterion tends to be met in the Northern Hemisphere, 
but not in the Southern Hemisphere, during midwinter. This difference 
may be related to the dissimilar midwinter breakdown climatologies. In 
the Southern Hemisphere, the structure varies between midwinter and 
spring so that the vortex tends to meet the instability criterion in 
spring when the vortex normally breaks down. 
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ABSTRACT OF DISSERTATION 

A COMPARISON OF THE STRUCTURE AND FLOW CHARACTERISTICS OF THE UPPER 
TROPOSPHERE AND STRATOSPHERE OF THE NORTHERN AND SOUTHERN HEMISPHERES 

The general circulations of the Northern and Southern Hemispheres 
are compared with regard to the upper troposphere and stratosphere 
using atmospheric structure obtained from satellite, multi-channel 
radiance data. Specifically, the data are from the Satellite Infrared 
Spectrometer (SIRS) instrument aboard the Nimbus 3 spacecraft. The 
inter-hemispheric comparisons are based on two months of data (one 
summer month and one winter month) in each hemisphere. 

The mean meridional circulation of the Southern Hemisphere strato- 
sphere in winter is studied using zonal ly-averaged vertical motion 
fields calculated by a heat budget technique. The mean pattern for 
July 1969 shows two main cells with the polar night jet stream associ- 
ated with an indirect cell. There are shown to be significant 
variations in the horizontal eddy heat flux and vertical motion patterns 
in the Southern Hemisphere in relation to the presence or absence of a 
minor midwinter warming. 

Horizontal eddy heat fluxes in the upper troposphere of both 
hemispheres are examined and in mid-latitudes found to be approximately 
equal in magnitude when averaged oveT the summer and winter month. The 
Southern Hemisphere mid- latitude eddy heat flux is also shown to have 
significant longitudinal variations, apparently associated with the 
location of the southern continents. 

The relative importance of standing and transient eddies in both 
hemispheres is investigated. The results confirm that, in general, 
standing eddies are less important in the Southern Hemisphere. This 
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distinction is true for the temperature structure as reflected in the 
eddy available potential energy (AE) , the eddy kinetic energy (KE) and 
the mid- latitude eddy heat flux. The distribution with latitude of the 
relative importance of standing and transient eddies is also studied. 

The zonal mean temperature and geopotential height gradients in 
the upper troposphere between 20° and 80° latitude are confirmed to be 
larger in the Southern Hemisphere in both the winter and summer, the 
inter-hemispheric difference being much greater in summer. Comparisons 
between the hemispheres of zonal and eddy available potential energy 
(AZ and AE) and zonal and eddy kinetic energy (KZ and KE) indicate that 
the ratios of AZ to AE and KZ to KE are larger in the Southern Hemi- 
sphere . 

In relation to the disparity in the breakdown climatologies of the 
Arctic and Antarctic stratospheric polar vortices, the structure of the 
two vortices is shown to be different so that a previously-derived 
instability criterion tends to be met in the Northern Hemisphere, but 
not in the Southern Hemisphere, during midwinter. In the Southern 
Hemisphere, the structure varies between midwinter and spring so that 
the vortex tends to meet the instability criterion in spring when the 
vortex normally breaks down. 


Robert Frederick Adler 
Department of Atmospheric Science 
Colorado State University 
Fort Collins, Colorado 80521 
April, 1974 
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CHAPTER 1 


PURPOSE AND STRUCTURE OF DISSERTATION 

The primary objective of this study is to investigate and compare 
particular aspects of the general circulations of the Northern and 
Southern Hemispheres, especially with regard to upper tropospheric and 
stratospheric layers. The study is based on atmospheric structure in 
both hemispheres derived from satellite, multi-channel radiance data. 
The satellite data provide information on atmospheric structure in both 
hemispheres with equal, or nearly equal, time and space resolutions, 
thus providing an excellent basis for inter-hemispheric comparisons. 

The major topics explored in this study include: 

a. atmospheric structure as reflected in the quantities zonal 
and eddy available potential energy, and zonal and eddy 
kinetic energy; 

b. upper tropospheric horizontal eddy heat flux and its spatial 
distribution; 

c. the relative importance of standing and transient eddies; 

d. the structure of, and the mean meridional circulation in, the 
winter stratospheric polar vortices. 

Previous observational studies of the Southern Hemisphere general 
circulation have been handicapped by a lack of horizontal resolution 
and by a relative lack of data in the upper troposphere and strato- 
sphere. Data from the International Geophysical Year (IGY) in 
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1957-1958 form the basis for most previous studies. Approximately 70 
radiosonde stations in the Southern Hemisphere formed the IGY upper-air 
network. This is essentially the same network that exists at present. 

In the Northern Hemisphere approximately 500 upper-air stations report 
regularly. Despite an effort to place radiosonde stations in data- 
sparse regions over Southern Hemisphere oceans, the IGY data set has 
very large gaps in its coverage. A post-analyses of the IGY data 
produced daily surface and 500 mb analyses which also have been used in 
a number of general circulation studies. Above 500 mb there are no 
post-analysis maps available and the quantity of upper-air data 
decreases with height, especially above 100 mb. However, some inves- 
tigations of the stratosphere in the Southern Hemisphere have been 
carried out using satellite radiance information to infer atmospheric 
structure. 

The present investigation is based on atmospheric structure 
derived from the Satellite Infrared Spectrometer (SIRS) instrument 
aboard the Nimbus 3 spacecraft. Temperature information is obtained 
using linear regression equations based on a comparison of satellite 
radiances and radiosonde data for the specific months and areas consid- 
ered in this study. Daily analyses are produced using twenty- four 
hours of satellite data. This daily data set has an adequate spatial 
distribution for hemispheric analysis of patterns of the sizes impor- 
tant in the large-scale energy and transport phenomena of the general 

circulation. Because of factors discussed in a following chapter, the 

i 

technique used to obtain atmospheric structure from the SIRS radiances 
is at its best in the upper troposphere and stratosphere. The satellite 
data, therefore, provide information concerning atmospheric structure 
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with adequate spatial resolution, with equal resolution in the two 
hemispheres and with the best results to be expected in the upper 
troposphere and stratosphere. Although there are limitations, the 
satellite radiance data have the potential to provide a great deal of 
information concerning the comparison of the circulations of the North- 
ern and Southern Hemispheres. To check the validity of computations 
made with the SIRS-based structure, a second set of calculations is 
carried out in the Northern Hemisphere using National Meteorological 
Center (NMC) gridded data. Although some comparisons between the 
hemispheres are made of quantities based on structure through the depth 
of the atmosphere, they are generally restricted to tropospheric layers 
above 500 mb, and the stratosphere. One summer month and one winter 
month are used in each hemisphere. 

For reasons discussed in Chapter 3 the structure obtained from the 
SIRS data tends to underestimate gradients. Mean meridional gradients 
are incorrectly estimated only slightly. However, variations in the 
east-west direction and the variance around a latitude circle of a 
variable, are significantly underestimated. Consequently, the satellite 
data would seem most appropriate for comparing mean meridional gradi- 
ents. However, these gradients can be, and have been, accurately 
estimated from monthly mean or even climatological data in a number of 
studies for both hemispheres. Terms which involve longitudinal 
gradients calculated on a daily basis are much more difficult to 
estimate from the limited radiosonde network of the Southern Hemisphere. 
Therefore, although there is a systematic underestimation of east-west 
gradients, the satellite data gives the time and space coverage so that 
a great deal of information about the eddy terms can be obtained. 
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Because of the underestimation, absolute magnitudes will not be as 
important as relative values in this study. 

Chapter 2 of this dissertation summarizes previous studies of the 
Southern Hemisphere's general circulation and investigations of inter- 
hemispheric comparisons. Chapter 3 describes the satellite system, the 
technique used to obtain atmospheric structure from the satellite 
radiance data and the limitations of the method. 

Chapter 4 contains the results of the inter-hemispheric comparisons 
in the troposphere. Values are presented for the zonal and eddy compo- 
nents of available potential energy (A) and kinetic energy (K) , with 
the emphasis being on the contributions of upper tropospheric layers .to 
the vertically integrated total. The eddy flux of sensible heat in the 
upper troposphere is also discussed. Spatial variations and the 
relative contributions of standing and transient eddies are examined. 

Differences between the hemispheres with regard to the winter 
stratosphere are investigated in Chapter 5. The relation of vortex 
structure to breakdown characteristics for the two hemispheres is 
discussed and calculations of the mean meridional circulation in the 
polar stratosphere of the Northern and Southern Hemispheres are presen- 
ted. Variations in wave characteristics and in the mean meridional 
circulation in the Southern Hemisphere are examined in relation to the 
presence or absence of a minor midwinter warming. 

A summary of the study is given in Chapter 6. 



CHAPTER 2 


REVIEW OF LITERATURE ON SOUTHERN HEMISPHERE GENERAL CIRCULATION 
AND DIFFERENCES BETWEEN THE HEMISPHERES 

2.1 Introduction 

Previous general circulation studies of the Northern and Southern 
Hemispheres reveal significant differences between them. This is not 
unexpected since there are great contrasts in their physical features. 

A good description of the meteorologically pertinent physical charac- 
teristics of the Southern Hemisphere is given by Taljaard (1972a) . 
Basically, three main factors contribute to the dissimilarity. First, 
a much smaller proportion of the Southern Hemisphere is covered by 
land. This is especially true of the middle latitudes. Between 40°S 
and 60°S, less than 5 % of the area is covered by land; in the same 
latitudes of the Northern Hemisphere, the percentage is approximately 
50%. Therefore, while around latitude circles of the mid-latitudes of 
the Northern Hemisphere there is a large land-sea contrast, in southern 
mid-latitudes ocean exists almost exclusively. Near the equator (less 
than 15 ° latitude), the percentage of land is approximately the same 
in the two hemispheres. Only in polar regions (greater than about 70° 
latitude) does the percentage of land in the Southern Hemisphere become 
significantly greater than in the Northern Hemisphere. This contrast 
in percentages at high latitudes reflects, of course, the presence of 
the Antarctic continent in the south polar region. This large, domed 
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land mass at the pole is the second main factor in the difference 
between the hemispheres. The third is the orography. The only moun- 
tain range in the Southern Hemisphere, which compares to those of the 
Northern Hemisphere, is the Andes Range of South America which runs 
south from the equator along almost the entire length of the west coast 
of the continent. Its highest portions are north of approximately 3S°S, 
and the range as a whole, and especially its southern portion, is 
relatively narrow. South Africa is a plateau averaging about 1000-2000 
meters, while Australia is a low plateau averaging below 500 meters in 
elevation. As a result, the Southern Hemisphere, with the exception 
of the Andes Range, does not possess orographic barriers of the number 
and lateral extent of those of the Northern Hemisphere. 

2.2 Radiation budget 

Observations of the earth’s radiation budget from satellites and 
calculations of the budget from radiation models both reveal significant 
contrasts between the Northern and Southern Hemispheres. However, 
budget parameters averaged over each hemisphere and over the year show 
no obvious differences between the hemispheres. Vonder Haar and Suomi 
(1971) summarize satellite measurements showing that the Northern and 
Southern Hemispheres have almost the same annual average albedo (30%) 
and emitted infrared radiation. But, while on a season-by-season basis 
the emitted infrared radiation values compare favorably in the two 
hemispheres, the albedo values do not. High albedos occur in the 
Northern Hemisphere in winter and spring and in the Southern Hemisphere 
in summer and fall. A very low albedo (23%) is observed for the 
Southern Hemisphere winter. This low seasonal value is related to the 
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high percentage of ocean coverage of the sunlit portion of the hemi- 
sphere during the polar night. 

More important to the general circulation than hemispheric values 
of radiation budget components are the latitudinal variations of these 
parameters. Vonder Haar and Suomi (1971), Vonder Haar (1972) and 
Raschke, Vonder Haar, Bandeen and Pasternak (1973) discuss pole-to-pole 
latitudinal profiles of radiation budget parameters as observed from 
satellites, Newell, Vincent, Dopplick, Ferruzza and Kidson (1970) 
compare satellite observations of net radiation as a function of lati- 
tude with various computations from radiation models on a seasonal 
basis. Sasamori, London and Hoyt (1972) compare their radiation model 
calculations for the Southern Hemisphere with satellite observations. 

In general, the observations compare favorably with the calculations, 
especially in regard to north-south gradients. 

The mean annual pole-to-pole profile of net radiation from 
satellite observations is given in Figure 2.1. On an annual basis, 
assuming no storage in the planetary system, the net radiation distri- 
bution is related to the northward transport of energy by the atmosphere 
and ocean by 

2 r4> i 

F(<h) = 2ira (4>) cos $ d$ 2.1 

' -ir/2 

where F(tf>i) is the northward transport or flux of energy across latitude 
i{>2, R^ is the net radiation, and a is the earth’s radius. The net 
radiation gradient is a sort of forcing function for the circulation of 
the atmosphere and ocean, although it is , because of its high dependency 
on clouds and thermal structure, not independent of the circulation. 
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Fig. 2.1 Mean annual, satellite-based net radiation distribution 
(from Vonder Haar, 1972) . 
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Fig. 2.2 Mean seasonal values of the meridional profiles of net 
radiation obtained from satellite measurements (from 
Vonder Haar and Suomi, 1971). 
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From Figure 2.1 it is obvious that the equator-to-pole gradient of 
R^ is larger in the Northern Hemisphere than in the Southern Hemisphere 
(Vonder Haar and Suomi, 1971). However, the equator-to-pole difference 
in net radiation may be misleading with regards to possible relations 
to the general circulation. That is, a stronger equator-to-pole net 
radiation difference does not necessarily mean that that hemisphere's 
general circulation is stronger in terms of mid- latitude energy fluxes. 
A better parameter perhaps is the hemispheric maximum of the required 
energy flux itself. This maximum will occur at the latitude of = 0. 
Calculations of the energy transport from net radiation considerations 
by Vonder Haar and Suomi (1971) and Newton (1972) indicate that the 
maximum in the two hemispheres is nearly identical, with the Southern 
Hemisphere having a slightly higher value. These similar values indi- 
cate that, on an annual balance, the circulation of the Southern 
Hemisphere must carry as much, or slightly more, energy poleward at 
the latitude of maximum transport. The energy flux is divided between 
the atmospheric circulation and the ocean circulation. Although 
calculations of the ocean contribution are uncertain. Sellers (1965) 
concludes that the percentage of the energy transport carried out by 
the oceans is greater in the Northern than in the Southern Hemisphere. 
From the estimates of total transport from the satellite net radiation 
measurements, it then follows that the atmosphere of the Southern 
Hemisphere must transport more energy polewards than its counterpart 
north of the equator. This conclusion is heavily dependent upon the 
evaluations of the ocean transport in the two hemispheres. As shown by 
Vonder Haar (1972) and Vonder Haar and Oort (1973), there are signifi- 
cant discrepancies between the ocean fluxes given by Sellers (1965) and 
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the ocean contribution calculated as a residual of measured atmospheric 
transport and satellite-measured net radiation. 

The annually-averaged storage term is usually neglected, even 
though there is a lack of good estimates of it. On a seasonal basis, 
however, the storage is probably larger and not negligible. Therefore, 
season-by-season, the net radiation is not balanced by the divergence 
of energy flux. However, the seasonal net radiation gradients (Figure 
2 . 2 ) are still important as a forcing for the general circulation. 
During the winter in both hemispheres, there is a reversal of the 
gradient near the poles. The reversal is more pronounced in the South- 
ern Hemisphere and shows up in the annual mean curve as the flat part 
of the curve near the South Pole in Figure 2.1. The net radiation 
gradients between high and low latitudes for winter are nearly the same 
for the two hemispheres. There is, however, a large variation between 
the hemispheres in the summer mean gradient. The equator-to-pole 
difference is only slightly larger in the Southern Hemisphere; however, 
it appears that the equator-to-pole gradient underestimates the hemi- 
spheric dissimilarity. While the value of the net radiation is nearly 
equal at the equator, in the latitude band 10°-40° there are substan- 
tially higher values of in the Southern Hemisphere. These higher 
values in the low latitudes, coupled with the lower values at higher 
latitudes (>70°) in the Southern Hemisphere, indicate a much greater 
latitudinal forcing in the Southern Hemisphere summer than in the 
Northern Hemisphere summer. It is also interesting to note that there 
is a major difference in the radiation budgets of the two polar regions 
in summer. While there is a negative net radiation poleward of about 
70°S, in the Northern Hemisphere the net radiation is positive at all 



latitudes, except zero at the pole. This disparity results largely 
from an inequality in albedo between the two regions. 

Longitudinal variations of components of the radiation budget may 
also be important to the atmospheric circulation. The geographical 
distributions of various budget components and the net radiation as 
observed from satellites are discussed by Vonder Haar and Suomi (1971) 
and Vonder Haar (1972). The mean annual pattern for the net radiation 
shows pronounced longitudinal fluctuations in the tropics, and a nearly 
zonal pattern in middle and high latitudes of each hemisphere. A 
comparison of the patterns at 15°N and 15°S shows longitudinal features 
of larger amplitude in the Southern Hemisphere (Vonder Haar and Suomi, 
1971). No such hemispheric difference can be detected at higher lati- 
tudes. Vonder Haar (1972) displays figures which seem to indicate 
that the Northern Hemisphere has a greater inter-annual variability in 
the longitudinal distribution of albedo and longwave radiation than the 
Southern Hemisphere. 

2.3 Mean tropospheric temperature structure and flow fields 

The recent publication of the Southern Hemisphere, upper-air 
climatology by Taljaard, van Loon, Crutcher and Jenne (1969) enables a 
comparison to be made for various levels with a Northern Hemisphere 
climatology such as that of Crutcher and Meserve (1970). The tempera- 
ture, height and wind climatologies of the Southern Hemisphere are 
discussed in detail by van Loon (1972a, b, c) and zonal mean cross- 
sections from pole-to-pole are discussed by Lorenz (1967), Newell et al. 
(1970) and van Loon (1972c). 
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The cross-sections of temperature and zonal wind show significant 
variations between the hemispheres. On a mean annual basis, the pel©° 
to-equator temperature gradient at 500 mb is larger in the Southern 
Hemisphere. This is also the case for the summer and winter seasons 
individually with the distinction between the hemispheres being small 
in winter and large in summer. This seasonal variation in hemispheric 
differences is also prevalent in the latitudinal distribution of net 
radiation (see Section 2.2). The distributions of net radiation and 
temperature are, of course, interdependent. 

Zonal mean geostrophic wind is related to latitudinal temperature 
gradient by the thermal wind relation. It follows that the maximum in 
the mean westerlies should be larger in the Southern Hemisphere on an 
annual average, with the difference much larger in the summer. The 
summer contrast is striking, with the Southern Hemisphere westerlies 
having a mean maximum speed about 50% higher than that found north of 
the equator. The difference in winter is less, with maximum speeds 
almost the same. 

Zonal cross-sections, however, do not reveal some of the most 
important inter-hemispheric dissimilarities. CrutcheT and Meserve 
(1970) are the most recent source for Northern Hemisphere, upper-air 
climatic maps, Comparison of the Northern Hemisphere 500 mb and 200 mb 
mean height fields of January and July with the Southern Hemisphere 
maps of Taljaard et al. (1969) for July and January respectively, reveal 
sharp differences with regard to zonal variations. The mean patterns 
in both seasons are much more zonal in the Southern Hemisphere, In the 
Northern Hemisphere a three-wave pattern is evident, especially in 
winter. The lack of large standing waves apparently reflects the lack 
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of geographic and orographic forcing in the Southern Hemisphere. 

However, the more zonal Southern Hemisphere patterns are not altogether 
symmetric about the pole. The distributions of SOO mb and 200 mb 
heights at 60°S for both January and July indicate that the polar 
vortex is not centered on the pole but offset towards about 50°E, 
towards Madagascar. This placement is related to a mean jet position 
at about 45°-50°S, centered at about 50°E, as evident on the geostroph- 
ic wind climatologies of van Loon, Taljaard, Jenne and Crutcher (1971). 
This is the major jet stream in the summer (January). In winter (July) 
the mean wind maps for 500 mb and 200 mb show two major wind maxima. 

One is at the aforementioned position at about 50°E. The other location 
at 500 mb is to the east of Australia at 180°E, at latitude 26°S. At 
200 mb the wind maximum has lengthened and cuts across the eastern 
part of the Indian Ocean, across southern Australia at about 30° S, and 
across the western Pacific Ocean. The highest analyzed isotach is for 
50 m sec” 1 . In summer (January) a weak maximum is also evident at 
200 mb across eastern Australia and the western Pacific Ocean at about 
30°S. The same feature is not distinctly evident on the mean map for 
500 mb. Therefore, while the Southern Hemisphere mean tropospheric 
flow is more zonal than that in the Northern Hemisphere, there are 
important zonal variations. 

2.4 Energy and momentum budgets 

The energy cycle of the Southern Hemisphere, and the components 
thereof, have not been studied in as much detail as those of the North- 
ern Hemisphere. Descriptions of the atmospheric energy cycle have 
therefore been based mainly on Northern Hemisphere studies such as 
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Oort (1964), However, some studies with a global scope have been 
conducted. Determination of the poleward transport of energy by the 
atmosphere and ocean combined on an annual basis has been accomplished 
using satellite net radiation data (see Section 2.2) or model computa- 
tions of the radiative budget of the planetary system. Some flux 
computations have been carried out in the Southern Hemisphere using 
IGY data. 

Newton (1972) gives a comprehensive discussion of the Southern 
Hemisphere general circulation in relation to global energy and momentum 
budgets. A consistent description of atmospheric and oceanic fluxes is 
arrived at by using absolute angular momentum, water and energy balance 
requirements along with flux computations where available. The momentum 
budget results discussed by Newton (from Newton, 1971) for both the 
annual case and the individual seasons indicate a number of differences 
between the two hemispheres. The Southern Hemisphere's tropics and 
mid- latitudes are, respectively, a greater surface source and sink of 
momentum for the atmosphere than the same areas of the Northern Hemi- 
sphere on an annual basis. The Southern Hemisphere also has, on an 
annual average, a greater poleward flux of momentum across 30 latitude. 
In addition, the momentum flux undergoes a much weaker annual variation 
in the Southern Hemisphere than in the Northern Hemisphere. However, 
maximum fluxes occur during winter both above and below the equator. 

Newton's discussion of the annual energy balance is based on the 
latitudinal distribution of net radiation values. These consist of a 
combination of calculated and satellite-observed values. The total en- 
ergy flux is divided into the latent heat flux, and the oceanic and 
atmospheric fluxes of potential heat which is defined as the sum of the 
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enthalpy and the potential energy. The latent heat fluxes were 
determined by precipitation-evaporation budget calculations. The flux 
curves determined by Newton are given in Figure 2.3. The oceanic 
component is based on an adjustment of the flux curve of Sellers (1965). 
Newton's modification decreases Sellers' estimate of the poleward 
transport by the ocean in the Southern Hemisphere, and increases the 
poleward transport by the ocean in the Northern Hemisphere. With the 
total energy transport (both the ocean and atmosphere) being nearly 
equal in the two hemispheres, the energy transport by the Southern 
Hemisphere atmospheric circulation must be greater than that of its 
northern counterpart . 

In the tropics (0-30°) the poleward potential heat transport by 
the atmosphere is much greater south than north of the equator. The 
latent heat flux in the southern tropics is also greater, but is 
directed equatorward in both hemispheres. In mid- latitudes the poten- 
tial heat flux and the latent heat flux are nearly equal in the two 
hemispheres. A breakdown of the atmospheric potential heat flux into 
toroidal and eddy components shows that in the tropics the variation 
between the hemispheres shown in Figure 2.3 is in the toroidal compo- 
nent. Thus, on an annual basis, the Southern Hemisphere has a much 
stronger Hadley cell in low latitudes. This stronger mean meridional 
circulation also accounts for the difference between the tropics of 
the two hemispheres with regard to the equatorward latent heat flux. 

In mid- latitudes the eddy component of the energy transport dominates 
and is approximately of the same magnitude in both hemispheres, accord- 
ing to Newton. 
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Fig, 2.3 Annual mean meridional flux of total energy for the 

earth-atmosphere system, and its apportionment between 
oceanic flux and atmospheric flux of potential heat and 
latent heat (from Newton, 1972), 
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Newell et al. (1970) also discuss the global energy budget of the 
atmosphere, and make global and hemispheric estimates of the components 
of the energy cycle in the troposphere (100-1000 mb) for the four 
seasons. The four basic components to the atmospheric energy cycle, 
zonal and eddy available potential energy (AZ and AE) and zonal and 
eddy kinetic energy (KZ and KE) , are calculated where possible. In 
addition, conversion rates, generation rates and dissipation rates are 
calculated. The AE and KE computations are severely hindered by the 
lack of data giving longitudinal variations on a day-to-day basis. The 
values for AZ and KZ indicate that for respective seasons, the Southern 
Hemisphere always has the larger values, the greatest difference 
between the hemispheres being during the summer. The quantities AZ and 
KZ are computed from climatological cross-sections. The energy con- 
version from AZ to AE is calculated using climatological values of 
temperature. During the winter months in each hemisphere, the AZ to AE 
conversion rates are of the same approximate magnitude; during the 
summer, the Southern Hemisphere conversion rate is much larger. The 
AE to KE conversion rate is not estimated in the Newell et al. study. 

2.5 Standing and transient eddies 

One of the major observed differences between the general circula- 
tions of the two hemispheres is the relative importance of standing and 
transient eddy motions. In the Northern Hemisphere, standing waves are 
evident in upper-air climatologies with mean troughs east of the Asian 
and North American continents and a weaker third trough in eastern 
Europe. No such strong stationary wave pattern is evident in the 
Southern Hemisphere climatologies, although the flow is not completely 
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symmetric about the pole (see Section 2.3). Van Loon and Jenne (1972) 
do show the presence of standing waves 1, 2, and 3 in the Southern 
Hemisphere, but the amplitudes of wave numbers 2 and 3 are much less 
than the same ones in the Northern Hemisphere. 

The relative lack of a standing eddy component of motion in the 
Southern Hemisphere is. also evident in calculations of heat and momentum 
transport and in calculations of eddy kinetic energy. Obasi (1963) and 
Kao, Jenne and Sagendorf (1971a) indicate that the transient eddies 
account for almost all of the poleward transport of momentum in the 
Southern Hemisphere. In addition, Obasi (1965) points out that what 
standing eddies there are in the Southern Hemisphere extract energy 
from the mean flow, unlike the standing eddies in the Northern Hemi- 
sphere which supply energy to the mean flow. Analyses of meridional 
eddy heat transport by Robinson (1970), Newton (1972) and Kao, Jenne 
and Sagendorf (1971b) in the Southern Hemisphere again reveal the 
overwhelming importance of the transient waves. On the other hand, 

Oort (1971) states that, in the Northern Hemisphere, the standing eddies 
transport half of the sensible heat during the winter months. 

Kao, Jenne and Sagendorf (1970) investigate the spectral properties 
of Southern Hemisphere eddy kinetic energy at 500 mb. Wave number 
spectra for both the total motion and for the stationary waves are 
given. A similar study (Kao and Wendell, 1970) gives the Northern 
Hemisphere version. For both the zonal velocity and the meridional 
velocity, the standing waves contribute less to the total eddy kinetic 
energy in the Southern Hemisphere than they do in the Northern Hemi- 
sphere. This difference is larger for the meridional velocity. 
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A comparison of spectra of kinetic energy for the Northern and 
Southern Hemispheres reveals another interesting hemispheric variation. 
The magnitudes of the zonal and meridional components relative to each 
other are different in the two hemispheres. Wooldridge and Reiter 
(1970) , using relative separation velocities obtained from GHOST balloon 
pairs, indicate a much larger ratio of meridional to zonal spectral 
density at synoptic-scale wavelengths in the Southern Hemisphere than 
is found in the Northern Hemisphere for similar lagrangian studies. A 
comparison of the wave number spectra of Kao and Wendell (1970) for the 
Northern Hemisphere and Kao, Jenne and Sagendorf (1970) for the Southern 
Hemisphere also indicates a dissimilarity in the spectral densities of 
zonal and meridional flow. In the Northern Hemisphere, the spectral 
density of u is larger than that of v at the synoptic scale (~ wave 
number 6), while in the Southern Hemisphere, the spectral density of 
v is larger. 

Although the studies just discussed emphasize the relative lack of 
standing eddies in the Southern Hemisphere, there are features in the 
cloudiness and precipitation climatologies (see van Loon, 1972d) which 
indicate some sort of geographic anchoring. In the western South 
Pacific Ocean, there is a narrow band of mean cloudiness and mean pre- 
cipitation running from the northwest to the southeast from a point at 
about 10°S, 1S0°W to 35°S, 1S0°W. The band is evident in both winter 
and summer, but is more obvious in summer. A second cloud band, 
similarly oriented, occurs east of South America; however, this band is 
not as distinct as the Pacific Ocean one. No similar cloud area exists 


in the Indian Ocean. 
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The reason for the existence and location of these bands is uncer- 
tain. The cloudy region in the Pacific Ocean is coincident with a 
maximum in cyclone frequency (Taljaard, 1972b) .^-Although the equator- 
ward limit of the clouds is associated-with the inter-tropical 
convergence zone and tropical cyclones, farther south the band is 
associated with frontal troughs and extratropical low pressure centers. 
In the Atlantic Ocean, the cyclone frequency maximum parallels, but is 
just south of, the cloud belt. The cause of this preference in cyclone 
tracks is unclear. 

The mean upper-level wind fields (van Loon, 1972c) show middle and 
upper troposphere wind speed maxima in the region of the Pacific Ocean 
cloud band. The axis of the jet stream is located at about 28°S. In 
January the maximum is delineated best at 200 mb as an extension east- 
ward from Australia of relatively high wind speeds with a peak of about 
23 m sec' 1 . In these longitudes there are actually two peaks, the other 
being farther poleward at about 50° S. In July at 200 mb there is a 
large jet stream center at 28°S centered over Australia and extending 
eastward with the speeds decreasing eastward. At 500 mb a maximum in 
the westerlies is located east of Australia at 28 S, 180 W. Therefore, 
in both seasons, there are upper-level features of importance occurring 
in association with the Pacific Ocean cloud band. 

The upper-level support for the Atlantic Ocean cloud band is 
indistinct, possibly because of a lack of upper-air data over the 
region. There are, however, indications on some of the isotach analyses 
that the band is located in the left-rear quadrant of a jet stream 
isotach maximum, which, in the Southern Hemisphere, should show upper- 
level divergence through vorticity considerations. 
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In the summer (January) , the Pacific Ocean cloud band extends 
southeastward from the inter-tropical convergence zone (ITCZ). Across 
the Indian Ocean and to the northeast of Australia to lSO^W, the ITCZ 
is located at approximately 10°S. East of about 140°W eastward to the 
Indian Ocean, the mean position of the ITCZ in summer is on the equator. 
Consequently, at the northern end of the cloud band, there is low-level 
convergence at the ITCZ. The southerly location of the ITCZ across the 
Indian Ocean and the western Pacific Ocean and the existence of a mean 
wind speed maximum in the middle and upper troposphere probably provide 
the conditions for cyclone development in the region of the cloud band. 

The extension of the ITCZ southward across the eastern half of the 
hemisphere in January is related to the strong northerly, low-level 
flow coming across the equator from the Asian winter monsoon. The 
upper- level flow in the Southern Hemisphere may also be affected by 
monsoon conditions on the other side of the equator. Therefore, in 
trying to find complete explanations for Southern Hemisphere climato- 
logical features, one must consider the global atmosphere as a whole. 

2.6 Stratosphere 

The stratospheres of the Northern and Southern Hemispheres are 
quite alike during the summer months and quite dissimilar during the 
winter. The circulation during summer in both hemispheres is dominated 
by a warm anticyclone centered over the polar regions. In the Southern 
Hemisphere, the center of the anticyclone is slightly off the pole 
towards 1S0°W (Labitzke and van Loon, 1972), The isotherms are nearly 
concentric around the pole. The Northern Hemisphere in July shows 
basically the same type of height and temperature patterns with the 



22 


anticyclone centered very close to the pole (Craig, 1965). The temper- 
atures in the Southern Hemisphere at 30 mb are higher than those in the 
Northern Hemisphere poleward of about 30°, and colder equatorward of 
that latitude, indicating a stronger meridional temperature gradient in 
the Southern Hemisphere (Labitzke and van Loon, 1972). The 30 mb 
heights in the polar region are, however, lower in the Southern Hemi- 
sphere because of a much colder troposphere. At 50 mb the height 
gradients indicate that stronger mean easterlies are present in the 
Northern Hemisphere. However, because the meridional temperature 
gradient is larger in the Southern Hemisphere, the easterlies increase 
with height faster there. Although Muench (1968) shows evidence of 
large-scale disturbances in the Northern Hemisphere summer stratosphere, 
both hemispheres are considered synoptically inactive at these heights 
relative to occurrences during the winter. 

From summer to winter in both hemispheres, temperature gradients 
in high latitudes reverse, westerlies take the place of easterlies, and 
the polar-night vortex is formed. Climatological ly, the Southern 
Hemisphere vortex has lower temperatures and lower heights. The south- 
ern vortex is also more symetric about the pole than its Northern 
Hemisphere counterpart. The climatic 50 mb height and temperature 
fields for the Northern Hemisphere (Craig, 1965) show pronounced wave 
numbers 2 and 3, while the flow at that level in the Southern Hemisphere 
is almost completely zonal at high latitudes. 

The major reason for the difference in the winter stratospheric 
climatologies is that the wintertime polar vortices of the Northern and 
Southern Hemispheres exhibit contrasting breakdown climatologies. In 
most years, the Northern Hemisphere vortex breaks down completely in 
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midwinter. Often the voTtex never becomes re-established. A detailed 
summary of events during Northern Hemisphere breakdowns is given by 
Reiter (1969a), The Antarctic vortex, on the other hand, has not been 
observed to undergo the complete, midwinter breakdown typical of the 
Arctic vortex, although midwinter minor warmings have been observed 
(Labitzke and van Loon, 1972). The final, complete collapse of the 
Antarctic vortex does not occur until spring. 

This difference in breakdown characteristics between the two 
hemispheres has not been completely explained. Case studies in the 
Northern Hemisphere have shown that the main source of energy for the 
breakdowns is an eddy flux of energy from the troposphere (Julian and 
Labitzke, 1965 and Matsuno, 1971). The waves which are important in 
the vertical flux into the stratosphere are the very long, slowly 
moving or stationary waves. In the Northern Hemisphere during the 
winter, stationary wave numbers 2 and 3 are prominent in the climato- 
logical fields of the upper troposphere. These stationary waves are 
associated with the distribution of geography and orography in the 
Northern Hemisphere. In the Southern Hemisphere mid- latitudes, the 
land-sea variation is essentially nonexistent, and the orographic 
barriers are minor compared to those in the Northern Hemisphere. The 
effect of this rather uniform surface is apparently a relative lack 
of standing eddies which may result in a lack of upward propagation of 
energy from the troposphere to the stratosphere in the Southern Hemi- 
sphere. Since this vertical flux of energy has been identified as the 
major source of energy for the Northern Hemisphere breakdowns, the lack 
of standing waves in the Southern Hemisphere may be related to the lack 
of midwinter breakdowns in the Southern Hemisphere. 
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A number of studies (Reed, Wolfe and Nishimoto, 1963} Mahlman, 

1966} Julian and Labitzke, 1965; Vincent, 1968) have shown that the mean 
meridional circulation in the lower half of the stratosphere in the 
Northern Hemisphere winter has a two-celled structure with rising motion 
in low latitudes and over the pole and sinking motion in mid- latitudes. 
This zonal ly-averaged vertical motion pattern exists before, during and 
after midwinter breakdowns of the polar vortex (Mahlraan, 1966), In the 
Southern Hemisphere the mean meridional circulation is practically 
unknown. Rubin and Weyant (1963) do indicate a sinking motion over the 
pole below 75 mb. Also Reiter (1969a, 1971), on the basis of ozone 
measurements over Antarctica, postulates a general sinking motion over 
the pole in the stratosphere. This descent over the polar region would 
tend to compensate for decreasing temperature due to the large radiative 
loss during the polar night. 

Some previous studies of Southern Hemisphere stratospheric circu- 
lation have used atmospheric structure inferred from satellite radiances. 
Nordberg, McCulloch, Foshee and Bandeen (1966), Julian (1967) and Shen, 
Nicholas and Belmont (1968) use 15 micron radiation data to interpret 
Southern Hemisphere stratospheric temperature patterns. Fritz (1970) 
and Fritz and Soules (1970) examine stratospheric temperature variations 
in the Southern Hemisphere using 15 micron radiances from the SIRS 
experiment. Smith, Woolf and Jacob (1970) and Miller, Finger and 
Gelman (1970) use 30 mb temperature and/or height information derived 
from SIRS data to study wintertime stratospheric patterns in the 
Southern Hemisphere, 



CHAPTER 3 


METHOD USED TO OBTAIN ATMOSPHERIC STRUCTURE FROM 
NIMBUS 3 SIRS DATA 


3.1 Satellite system 

The primary data used in this study are from the Satellite 
Infrared Spectrometer (SIRS) experiment aboard the Nimbus 3 spacecraft. 
The satellite orbit is essentially circular at 600 nautical miles above 
the surface, is sun synchronous with a local noon equator crossing, and 
has an inclination of 100° which results in the orbital track never 
being poleward of latitude 80°. Successive orbits are 26° of longitude 
apart at the equator, and the period of the orbit is 107 minutes. 

The SIRS instrument measures the infrared radiation from the earth 
in eight spectral bandpasses, ranging from 11 to 15 microns (899 to 
669 cm -1 ). The eight channels are thus in a spectral interval where 
the atmospheric absorption varies from near 0% in the atmospheric 
"window" at 11 microns to near 100% at 15 microns in the CO 2 absorption 
band. This change in absorption with wavelength is the basis for 
obtaining temperature information as a function of height from the SIRS 
radiances. The field of view of the instrument is 11.5° on a side, 
which results in a view of the surface of 120 nautical miles on a side. 
A detailed description of the SIRS instrument is given by Wark, 
Hilleary, Lienesch and Clark (1969). 
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A number of techniques to obtain atmospheric thermal structure 
from multi-channel, satellite radiance data have been developed and 
used. They are reviewed by Fleming and Smith (1972). The methods can 
be separated into two categories: matrix inversion and iteration. 

Although some require no prior knowledge of atmospheric thermal struc- 
ture, most are based on additional information given by atmospheric 
temperature and radiance statistics. The simplest statistical -based 
scheme is a regression technique, in which near simultaneous satellite 
radiance observations and radiosonde observations are compared to 
develop regression equations, with atmospheric temperature as the 
dependent variable, and radiances as the independent variables. This 
type of procedure applied to Nimbus 3 data is discussed by Smith, Woolf, 
and Jacob (1970). The method used to obtain the thermal structure for 
the present study is similar to the regression technique of Smith, 

Woolf and Jacob (1970), but uses the geopotential height difference 
between pressure levels as the dependent variable. 

3.2 Physical basis for thickness - radiance regression technique 
The thickness D between pressure surfaces in the atmosphere is 
proportional to a weighted mean temperature in the layer, as given by 
equation 3.1. 

R f Pl 

D=z 2 -z i =— T dlnp 3.1 

g J P2 

The radiance observed by a satellite in a particular channel is also a 

function of a weighted mean temperature through a layer. Equation 3.2 

is the infrared radiative transfer equation. is the radiance 
til 

observed in the i channel by the satellite instrument, W^ is the 
Planck function defined in equation 3.3, is the i^ channel 
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transmittance function, and the subscript s represents surface condi- 
tions . 

.Pc 


, rs d-r. 1 

N. = - IV. -rs——- dlnp + - W. (T ) 
i 7 r J i dlnp r tt i s 


f P S dT, 


J 


dlnp 


i- dlnp 


3.2 


c l k i 

“i * e tc 2 k i /T) - 1 


3.3 


dt. 

l . 


In equation 3.3, ^ is the wave number. The term t q— is the weighting 
function for the i^ channel. The weighting functions for the eight 
SIRS channels are shown in Figure 3.1. Channels near the middle of the 
15 micron absorption band have maxima at high altitudes, while channels 
out on the wing of the absorption band have peaks at lower altitudes. 
Channel 1 (899cm' 1 ) is in the atmospheric "window." 

For the moment, assume a cloudless atmosphere and neglect the 
second term on the righthand side of equation 3.2, the surface contri- 
bution term. This term would be identically zero when dealing with 
channels 4 through 8 where the surface elevation is small. The term 
can also be considered known because T^ could be calculated using 
channel 1 (899cm' 1 ), and W^(T g ) could then be computed using equation 
3.3. With the surface contribution term neglected, equation 3.2 
becomes 


1 r P s dt . 

N i'7j o *i3IH? dln P 3.4 

The Planck function kh can be closely approximated by a linear function 
of temperature over a limited range of temperature when dealing with 
the 11-15 micron portion of the spectrum. For example, over a 40°K 
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Fig. 3.1 Derivative of transmittance with respect to the logarithm 
of pressure. These functions approximately describe the 
relative sensitivity of the eight SIRS radiance observa- 
tions to temperature variations in various altitude layers 
of the atmosphere (from Smith, Woolf and Jacob, 1970). 
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temperature range from 214 I> K to 254°K, a linear fit would reproduce 
Planck function values within 3% for the wavelength range from 13.4 
micron (750cm” 1 ) to. 15.0 micron (670cm” 1 ). Therefore, 

W^= C 3 + Ci+T 3.5 


Substituting equation 3.5 into equation 3.4, one obtains 


. 1 


N. = - c 3 
1 n 3 


s dx 
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dx. 

1 

dlnp 


dlnp 


3.6 


The first term on the right side of equation 3.6 is known, and sub- 
tracting it from hh results in 


n : = - c 4 

1 IT n 


,P S dx . 

T 3ISF dl "» 


3.7 


where N.’ = N. - — Co 
1 1 tt 0 


s dx . 


dx , 


dlnp. If were a rectangular function 


as shown in Figure 3.2, equation 3.7 would reduce to 


P2 

NJ = y c 4 T dlnp 3.8 

Jp i 

In other words, the N? would be proportional to a weighted, layer-mean 
temperature between two pressure levels, which in turn is proportional 
to the thickness D of the layer. Unfortunately, a glance at Figure 3.1 
shows that the weighting functions are not rectangular. However, it 
is also evident from an inspection of Figure 3.1 that a linear combina- 
tion of the weighting functions could be made that is "more rectangular" 
than any individual weighting function. Thus, a linear combination of 
radiances would be more closely proportional to the layer mean 




di*p 


Fig. 3.2 Idealized weighting function. 
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temperature or thickness than an individual radiance. Fleming (1972) 
discusses a mathematical technique to determine the optimum linear 
combination of radiances. In the present study, the linear combination 
is determined through a step-wise regression technique. 

3.3 Thickness regression method 

From the discussion in the previous section, one would expect a 
high correlation between the thickness D and a linear combination of 
channel radiances. To determine the coefficients of the linear combina- 
tion, a step-wise, least-squared error, multiple regression technique 
is used. The regression equation is 

D = ci0 + a^Ni + 012 N 2 + +agNg , 3.9 

where D is the thickness and through Ng are the satellite radiances 
in the eight channels. The dependent variables in the regression pro- 
cedure are the thicknesses for the following layers: 100-700 mb, 

700-500 mb, 500-300 mb, 300-200 mb, 200-100 mb, 100-50 mb, 50-30 mb, and 
30-10 mb. The independent variables are the radiances of the eight SIRS 
channels. Separate regression equations are determined for each layer, 
for each of the latitude zones 20-40°, 40-60°, and 60-80° in both hemi- 
spheres and for each month studied. The division by latitude is required 
because the relationship between thickness and radiance derived in 
Section 3.2 is based on the Planck function being approximately linear 
with respect to temperature. This approximation is only valid over a 
limited temperature range; therefore, the data is divided into latitude 
zones to restrict the temperature variation in each data group. 

The equations in Section 3.2 are for a cloud-free atmosphere. The 
regression technique is based on a comparison of thickness information 
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and cloud- free radiance data. When the derived regression equations 
are applied to other radiance data, only cloud-free data are used. To 
eliminate cloud-contaminated data from the comparison data set and from 
the application of the regression equations, a simple, objective "cloud 
check" procedure based solely on the satellite radiance information is 
used. The SIRS channels with weighting function peaks in the low 
troposphere (channels 1 through 3) are the ones most affected by the 
presence of clouds. Channels with weighting function peaks in the 
stratosphere (channels 7 and 8) are only rarely affected by the presence 
of tropospheric clouds. The radiance in the atmospheric window 
(channel 1, k = 899 cm -1 ) in the absence of clouds is related to the 
surface temperature through equation 3.3. In the presence of an over- 
cast, the channel 1 radiance is a function of cloud- top temperature, 
which is lower than the surface temperature. Therefore, channel 1 
radiances much below normal would indicate the presence of clouds; 
however, the radiance at this wavelength is highly variable even in 
the absence of clouds because of its dependence on the surface tempera- 
ture. Large changes occur along satellite tracks in the presence of 
sharp changes in surface characteristics such as land-sea boundaries. 
Also, large diurnal changes are present in the window radiance because 
of large changes in land surface temperature. Channels 2 and 3 are 
also affected by clouds, but are not so severely affected by surface 
characteristics as channel 1. Therefore, they are used for the "cloud 
check." 

SIRS radiance data is eliminated as being cloud- contaminated when 
the observed radiance in both channels 2 and 3 is below critical magni- 
tudes. The criteria for both channels were determined by a comparison 
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of SIRS data and satellite video data and are a function of latitude 
and season. The latitude profiles of the critical values for each 
season are given in Figures 3.3 and 3.4, Any data with both channel 2 
and 3 radiances below the critical levels shown in Figures 3.3 and 3.4 
are eliminated as cloud-contaminated. The units of radiance along the 

ordinate in the figures are ergs cm -2 sec -1 steradian -1 (cm -1 ) * 

abbreviated as (ergs/) . The "cloud check" procedure is simple and 

certainly not foolproof. Some slightly cloud-contaminated data may 

still find its way into the final data set. The procedure outlined 

above removes, on the average, about 15% of the original data points. 

and 3 radiances below the critical levels shown in Figures 3.3 and 3.4 
are eliminated as cloud-contaminated. The units of radiance along the 
ordinate in the figures are ergs cm -2 sec -1 steradian' 1 (cm' 1 ) , 

abbreviated as (ergs/). The "cloud check” procedure is simple and 
certainly not foolproof. Some slightly cloud-contaminated data may 
still find its way into the final data set. The procedure outlined 
above removes, on the average, about 15% of the original data points. 

To obtain the regression coefficients, comparison sets of radiance 
and thickness data are developed. The observations are near-simulta- 
neous and at, or close to, the same geographic position. Somewhat 
different techniques are used for the Northern and Southern Hemispheres. 
In the Northern Hemisphere, for the layer 100 to 1000 mb, the SIRS 
radiance data is matched with the National Meteorological Center (NMC) 
gridded data. For satellite data occurring within three hours of NMC 
map time, the NMC grid point data are linearly y interpolated to the 
satellite track position. For the layers above 100 mb, station data are 
matched with the satellite radiance information. The station location 
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Fig. 3.3 The distribution with latitude of critical values of 

channels 2 and 3 for the cloud check procedure for July 





60 ° 80 * 
South 


40 ° 


itical values of 
procedure for 






36 


must be within 180 nautical miles of the satellite position, and the 
satellite time must be within three hours of the station observation 
time. 

In the Southern Hemisphere, a different technique is used. Station 
data are matched with satellite radiances for all levels. Because of 
the relative lack of stations in the Southern Hemisphere and the need 
for a sufficient number of data points for the regression analysis, the 
match up restrictions are somewhat relaxed from those in the Northern 
Hemisphere. For each 24-hour period surrounding the 0000 GMT and 1200 
GMT radiosonde times, radiance information is placed into a five degree 
latitude by five degree longitude grid for the latitudes 20 o S to 80°S. 
The radiance fields are analyzed and the radiances in each of the eight 
channels are linearly interpolated from the grid locations to the 
stations. The match-up is restricted to times and locations when there 
are original satellite data within a geographic box 600 nautical miles 
on a side surrounding the station. This method is used for all levels. 
In the Southern Hemisphere no regression equations are developed for 
the 10-30 mb layer, because of a lack of radiosonde data at that alti- 
tude. The Northern Hemisphere regression equations for this layer and 
for the respective season are used to complete the Southern Hemisphere 
temperature structure. 

The choice of the regression equations to be used in computing 
thicknesses is based on certain statistical and practical criteria. 

The step-wise, least-squared error, multiple regression technique uses 
an F test to determine if additional independent variables should be 
introduced into the particular regression under consideration. In 
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addition, the number of independent variables is restricted to a maxi- 
mum of four. Most of the regression equations have three independent 
variables; however, a few have less than three. 

3.4 Regression statistics 

The regression equation coefficients for thickness used in this 
study are given in the Appendix along with the percentage of variance 
explained and the standard error of estimate (approximately equal to 
root mean square error) for each of the regression equations. The 
standard errors for each of the eight layers are given in Table 3.1 for 
each hemisphere and season. These statistics are for the three latitude 
bands combined (area weighted), and for the dependent data on which the 
regression equations are based. The errors are largest in the lowest 
layers because not quite all cloud-contaminated radiance information is 
removed by the "cloud check" procedure. In order to be absolutely sure 
that no cloud-contaminated data are entering the system, one would have 
to eliminate a large portion of the satellite radiances, allowing only 
the small remaining portion to be available for computations. As a 
result the "cloud check" parameters are set so that most cloud-contami- 
nated data are omitted, but some borderline cases are accepted; hence 
the larger errors are in the lowest layer. 

In general, the standard errors of estimate in Table 3.1 are 
comparable to those in other studies using a regression method. Smith, 
Woolf and Jacob (1970) show errors ranging from 1.0°C to 2.6°C for a 
May-June period using data from the same SIRS instrument. These values 
are similar to those in the first column of Table 3.1. Ohring (1972) 
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N. H 
July 69 

N. H 

January 70 

S. H 
July 69 

S. H 

January 70 

10-30 mb 

1.3° C 

2.1° C 



30-50 mb 

1.4 

3.0 

1.7° C 

1.8° C 

50-100 mb 

1.6 

2.3 

1.7 

2.2 

100-200 mb 

1.9 

2.0 

1.9 

2.4 

200-300 mb 

1.9 

2.4 

2.3 

1.9 

300-500 mb 

2.0 

2.8 

2.3 

2.0 

500-700 mb 

2.0 

3.5 

2.8 

2.1 

700-1000 mb 

2.6 

5.2 

3.2 

3.1 

Table 3.1 

. Hemispheric standard errors 
regression equations 

of estimate for thickness 

obtains errors 

of between 2°C 

and 3°C using 

radiances from 

the Nimbus 3 


Infrared Interferometer Spectrometer (IRIS) and a step-wise regression 
technique. 

In the Northern Hemisphere, standard errors are generally larger 
in the winter month. This difference is associated with the tighter 
gradients and more rapid movement of systems in the winter. Smith, 
Woolf and Jacob (1970) also observed an increase in standard error from 
summer to fall. Because of the larger errors, the regression results 
for January 1970 are examined in more detail for the 100-1000 mb layer. 
The January data is divided into three subsets: 1-6 January, 14-19 

January, and 26-31 January, Regression equations are derived for the 
subset of 1-6 January and are applied to the two independent sets. 
Separate equations are used for each of the latitude bands discussed in 
the last section. Root mean square errors are given in Table 3.2. The 
errors for the independent data are larger than those for the dependent 
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group, but there is, in general, less of a change between the two 
independent sets. Thus, the regression equations appear stable. When 
the three sets of data are combined and new regression equations deter- 
mined, the errors associated with the combined data drop to near the 
1-6 January set values. The above was done with regression equations 
restricted to a maximum of three independent variables or channels. 





Dependent 
1-6 January 
1970 

Independent 
14-19 January 
1970 

Independent 
26-31 January 
1970 

3 sets 
combined 

100- 

200 

mb 

1.9° C 

2.2° C 

2.6° C 

2.1° C 

200- 

300 

mb 

2.4 

2.8 

2.8 

2.5 

04 

O 

o 

r 

500 

mb 

2.7 

3.0 

3.1 

2.8 

500- 

700 

mb 

3.8 

4.8 

5.1 

3.8 

700- 

1000 mb 

5.2 

7.2 

6.6 

5.7 


Table 3.2. Root mean square errors for regression coefficients 
stability test 

3.5 Limitations of method 

The satellite-derived thickness fields used in this study are 
computed by applying the regression equations to radiance data (eight 
channels) at 5° latitude by 5 ° longitude grid locations. The data are 
processed and placed into the latitude-longitude grid system in the 
following way. Each set of eight observations of radiance (there are 
approximately 3000 sets for each hemisphere between 20° and 80° latitude 
for each 24 hours) is first checked for blatant errors and then for the 
presence of clouds as discussed in Section 3.3. If there are obvious 
errors or the data are determined to be cloud-contaminated, the set of 
observations is excluded. These two checks eliminate about 20% of the 
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data. Each accepted set of radiances is positioned at the nearest 5° 
by 5° grid location. Multiple sets of data at a particular location 
are averaged. For each 24 hour period an average of 400 of the 936 
possible locations contain original data. The remainder of the loca- 
tions are filled in by linear interpolation with a greater weight 
placed on the east-west direction. Geopotential height information is 
obtained by adding the thicknesses to a 1000 mb height field. This 
results in height data at the following levels: 700 mb, 500 mb, 300 mb, 

200 mb, 100 mb, 50 mb, 30 mb, and 10 mb. In the Northern Hemisphere 
NMC 1000 mb height fields are used. In the Southern Hemisphere 1000 mb 
height analyses are obtained by converting daily surface pressure 
analyses of the Commonwealth Bureau of Meteorology, Australia. The 
Southern Hemisphere surface pressure analyses are nearly hemispheric in 
scope, but where they are blank the field is filled in subjectively by 
a combination of extrapolation and climatology. 

The thickness and height fields produced in this study from the 
SIRS data have certain limitations and possible biases. Because only 
"cloud-free 1 ' data are. used in the satellite-based analyses, there is a 
bias of clear areas over cloudy areas which could produce a bias in the 
thickness of height fields if cloudiness is correlated with temperature 
and/or height. This correlation is certainly true in the low layers 
where low pressure centers are accompanied by areas of cloudiness. In 
the upper troposphere however, the correlation is much less. Although 
the failure to determine atmospheric structure in cloudy areas means 
that atmospheric eddies are not defined perfectly, the data on either 
side of cloudy areas are usually sufficient for adequate analysis. 



41 


In the following chapters, comparisons are made of general circu- 
lation parameters calculated in the Northern and Southern Hemispheres 
from SIRS-based atmospheric structure. In the Northern Hemisphere, 
calculations made with SIRS-based structure and ones made with NMC 
fields are compared. These comparisons reveal a systematic difference 
between the two sets of computations. Gradients in the NMC fields are 
underestimated in the SIRS-based fields. The amount is, in general, 
small with regard to meridional gradients and substantial with regard 
to longitudinal gradients. For example, the mean temperature gradient 
between 20°N and 80°N in the troposphere tends to be underestimated by 
about 10%. On the other hand, the variance of temperature around a 
latitude circle in the upper troposphere is underestimated by about 
50%, indicating that the amplitude of eddies is evaluated incorrectly 
by about 30%, since the variance involves squares of deviations. 

The underestimation of meridional gradients is minimal because 
separate regression equations are used for each of the three latitude 
bands. The variance of temperature around latitude circles is under- 
estimated because of a number of factors. First, some features, usually 
small, are not detected by the satellite-based structure. This is a 
relatively minor effect and is closely related to the distribution of 
satellite data. More important is the general result that the SIRS- 
based thickness fields tend to be "smoother” than the NMC fields which 
is partially due to the way the usable radiance data is distributed in 
space for a 24-hour period. Although a particular feature may be de- 
tected by the satellite data, the lack of a pass directly over the 
center of the feature results in an underestimation of the amplitude. 
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In addition, 24 hours of information are combined to produce the daily 
fields from which the SIRS-based calculations are made. These calcula- 
tions are compared to conventional data fields (NMC) at a particular 
time (1200 GMT). If the NMC fields were first averaged over 24 hours, 
the calculated values of temperature variance would than be less because 
the averaging would reduce the amplitudes of moving waves. The regres- 
sion technique also tends to aid in producing smoother fields than the 
NMC fields. The thickness-radiance regression equations work very well 
near the mean thickness for that latitude band, but tend to underesti- 
mate deviations from the mean, producing slightly "smoothed-out" fields. 

In addition to the effects discussed above which apply at all 
levels, in the lowest layers the distribution of satellite data, the 
slight smoothing of the radiance data and the effect of eliminating 
cloud-contaminated data combine to produce temperature fields which do 
not indicate the large east-west gradients that are present. The two 
calculations of the temperature variance agree best at higher eleva- 
tions. This is partly due to the decreasing effect of possible cloud 
effects, but more importantly due to the change in the character of the 
dominant waves with height. In the troposphere smaller, transient waves 
contribute significantly to the temperature variance; in the strato- 
sphere larger, quasi-stationary, or slowly moving, waves dominate. The 
large waves in the stratosphere are much more easily detected and more 
accurately depicted because of the time and space distributions of the 
satellite data. 
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Because of the limitations of the SIRS-based structure described 
in this section, the emphasis of this study is not on the absolute 
magnitude of quantities in each hemisphere, but on relative magnitudes 
and on parameters well estimated from the satellite-derived data fields. 



CHAPTER 4 


COMPARISONS BETWEEN THE HEMISPHERES IN THE TROPOSPHERE 

4.1 Available potential and kinetic energies 

4.1.1 Formulae and notation 

The comparison between hemispheres in this section will be made 
within the context of calculations of available potential and kinetic 
energies. The concept of available potential energy (A) and its role 
in the atmosphere's energy cycle have been discussed by Lorenz (1955, 
1967), Dutton and Johnson (1967) and Reiter (1969). The formulations 
for A and its components, the zonal available potential energy (AZ) and 
eddy available potential energy (AE), are essentially the same among 
the various authors. The so-called exact formula shows A to be a 
function of the variance of pressure on isentropic surfaces. An 
approximate formula for A (or AZ and AE) involving the variance of 
temperature on constant pressure surfaces is often used in calculations 
of available potential energy and is also adopted in this study. 

The notation used in this study closely follows that of Reiter 
(1969a, 1969b). Brackets represent an average over the coordinate 
given by the subscript. Parentheses indicate deviations from the mean. 
For example. 


where X is longitude. 


X = [X], 
The term (X) 


+ (X) x 4,1 

is the deviation of X from the 



45 


zonal average, [X]^. An average over both latitude and longitude is 
defined as 


[X] 


X,<t> 1 - s in4> 


[XJ X cos<j> dp. 


4.2 


so that [X] is an area-weighted mean. In equation 4.2, $ is the 
latitude confining the integration. If the integration covers the 


globe, <f = - / 0 . The expression 


™A>* = [X1 X - [X] X,^ 4 - 3 
defines the deviation of a zonal ly-averaged parameter from an average 
over both latitude and longitude. For a more detailed description of 
this notation the reader is referred to Reiter (1969b), 

With this notation the computational formulas for AZ and AE, 
respectively, are 


AZ 


“PI,)}]*., 



dp. 


4.4 


AE 


fPs 


d P» 


4.5 


where p is the surface pressure. The term [a], , is the mean static 

S A , 9 

stability given by 


[o] 


A, 4 


- [t(-£- + — )i =[_&£, Hi 

1 l c 9z JJ A,(f) 1 R e 3p J A ,< 


4.6 


In the calculations presented here, the temperature T in equations 
4.4 and 4.5 is determined from the thickness between constant pressure 
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surfaces. Strictly speaking, the thickness is proportional to the 
virtual temperature, not just the temperature. However, in all calcu- 
lations made in this study, no distinction is made between the two. 

The integrations of equations 4.4 and 4.5 are approximated by a summing 
over layers in the vertical. The static stability [a] ^ ^ is obtained 
using calculated mean temperatures at the upper and lower edges of the 
particular layer. For example, fcr] for the 300-500 mb layer is 
calculated using hemispheric averaged 300 mb and 500 mb temperatures. 

In the satellite-based calculations, these temperatures are obtained 
using regression equations derived in addition to those for thickness. 
The calculations are contained area-wise by latitudes 20° and 80° in 
each hemisphere. Vertically the integration is from the surface to as 
high as 10 mb. 

The total kinetic energy (K) of the atmosphere is usually divided 
into two components, zonal kinetic energy (KZ) and eddy kinetic energy 
(KE) , the formulas for which are given in equations 4.7 and 4.8 respec- 
tively. 


KZ 



fK « wjLi 

2g 


dp 


4.7 


KE 


f p s [(«)» * (V)*] X- , 

51 dp 

■'o 


4.8 


The resulting units are in terms of energy per unit area. The geo- 
strophic approximation is used to calculate the winds above in the 
present study. The term [v]^ cannot be evaluated because it is 

identically zero under the geostrophic assumption. However, [v], is at 

A 
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least one order of magnitude less than [u]^. Therefore, [u]^ >> [v]^ 
and KZ can be accurately calculated knowing just the mean zonal flow. 

4.1.2 Mean meridional gradients 

The zonal available potential energy (AZ) is a function of the 
variance of zonal ly-averaged temperature about [TL ., the area- 
weighted temperature usually calculated on a global or hemispheric 
basis. AZ is therefore very dependent on the mean meridional temperature 
gradient. Figure 4.1 shows the monthly mean, layer-mean temperature in 
the 300-500 mb layer in both hemispheres for January and July. For the 
Northern Hemisphere (parts a and b) , three curves are shown: climatology 

(based on data from Oort and Rasmusson, 1971), and two curves based on 
NMC analyses and SIRS-based analyses respectively. For the Southern 
Hemisphere (parts c and d) , there is a climatological curve (based on 
data from Taljaard, van Loon, Crutcher and Jenne, 1969) and the SIRS- 
based curve. All temperatures are computed from thicknesses. The 
climate curves are based on climatological 300 mb and 500 mb heights. 

The diagrams of Figure 4.1 display the differences between the 
hemispheres already noted in Section 2.3. In both seasons shown, the 
mean temperature gradient between 20° and 80° latitude is larger in the 
Southern Hemisphere. The inter-hemispheric difference is largest in 
the summer. In winter (Figure 4.1a and c) the gradient from 20° to 45° 
latitude is about the same with the climatological curves actually 
showing a sharper drop-off with latitude in the Northern Hemisphere. 
Poleward of 45° however, there is a much larger mean gradient present 
in the Southern Hemisphere. In summer (Figure 4.1b and d) , in addition 
to the obvious variation in mean gradient, there is a very distinct 



300-500 mb 



Fig. 4.1 Mean meridional temperature gradients in the upper 

troposphere (300-500 mb) in the Northern and Southern 
Hemispheres for January and July. 
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Fig. 4.1 Continued. 
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difference in low latitudes (between 20° and 30° latitude). In the 
Northern Hemisphere, there is little or no decrease in temperature 
between 20°N and 30°N; in the Southern Hemisphere there is a sharp 
decrease. This variation is probably related to the much stronger 
summer monsoon regime in the Northern Hemisphere. 

The SIRS-based curves in Figure 4.1 in general agree well with the 
NMC-based curves (in the Northern Hemisphere) and climatology. One 
significant difference occurs in Figure 4.1a between S0°N and 60°N. The 
maximum deviation between the SIRS-based temperature and the NMC tempera- 
ture is 1.5°K and is located at 55°N, At other latitudes, the NMC and 
SIRS curves closely agree. A similar difference does not show up in 
the Northern Hemisphere July (Figure 4.1b) diagram. In addition, neither 
of the Southern Hemisphere SIRS curves show an abrupt change of tempera- 
ture gradient as is evident at 55°N in Figure 4.1a. The probable 
reason for the overestimation of the temperature at 55°N by the SIRS- 
based technique is related to the distribution of satellite data in 
January 1970 in the Northern Hemisphere. That month in the Northern 
Hemisphere has substantially (- 25%) fewer original satellite data 
points than the other three data samples examined. In addition, the 
distribution of data with respect to longitude is uneven, with a relative 
lack of data in two regions of longitude, the primary one being over the 
western Pacific Ocean. This results in a greater underestimation of the 
amplitude of a temperature trough than the amplitude of the accompanying 
temperature ridge. Therefore, an overestimation of the zonal ly-averaged 
temperature occurs. 

The zonal kinetic energy (KZ) is closely tied to zonal ly-averaged 
geopotential height distributions through the geostrophic wind relation. 
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Figure 4.2 shows various comparisons of the height distribution at 
200 mb. The comparison in Figure 4.2a of NMC and SIRS curves reveals 
a substantial difference between the two curves in the latitude zone 
50°-60°N, with the SIRS values of [z]^ being larger by up to 95 m at 
55°N. This overestimation of the geopotential height is related to the 
overestimation of tropospheric temperature at the same latitudes as 
discussed in the last paragraph. None of the SIRS curves in Figures 
4.2b-d indicate a similar effect. 

In general, the latitude distributions of [z]^ in Figure 4.2 
display the same inter-hemispheric differences and seasonal variations 
that were evident in the upper tropospheric temperature distributions 
seen in Figure 4.1. For the winter month (Figures 4.2a and c) , the main 
difference occurs poleward of 45° latitude. At 75° latitude, the 
Southern Hemisphere climatological height is 350 meters lower than that 
in the Northern Hemisphere. In summer the difference at 75° latitude 
increases to almost 600 meters. For both months in the Southern Hemi- 
sphere, the SIRS-based curves of 200 mb [z]^ lie very close to the 
climatological curves. This agreement indicates that, at least in 
regard to mean zonal features, the months used in this study for the 
Southern Hemisphere are not atypical. 

Figure 4.3 shows the mean geostrophic zonal flow at 200 mb as a 
function of latitude in the Southern Hemisphere. The SIRS and climate 
curves, calculated from the geopotential height distributions in Figure 
4.2, agree fairly well. The jet stream moves poleward from July to 
January and the peak wind speed is larger in July. Although at 500 mb 
the peak wind is higher in January (summer) than in July (winter) 
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Fig. 4.2 Mean meridional geopotential height gradients at 200 mb 
in the Northern and Southern Hemispheres for January and 
July. 



Fig, 4.2 Continued. 
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Fig. 4.3 Distribution with latitude of zonal mean geostrophic 

wind at 200 mb in the Southern Hemisphere for July and 
January. 
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(van Loon, Taljaard, Jenne and Crutcher, 1971), the opposite is true 
at higher levels such as 300 mb and 200 mb, as seen in Figure 4.3. 

4.1.3 Results of available potential and kinetic energy comparisons 

The results of the computation of AZ are given in the first two 
lines of Table 4.1. Calculations are made on a daily basis, then 
averaged over the month. The contributions to the total AZ from the 
100-1000 mb layer, AZ (100-1000 mb), and the 300-500 mb layer, 

AZ (300-500 mb) are shown. Calculations for the Northern Hemisphere are 
presented using both the NMC fields and the SIRS-derived structure. 

The results of the AZ computations indicate the same hemispheric 
differences that were indicated in Figure 4.1. The summer disparity is 
very striking, with the Southern Hemisphere having twice the 
AZ(100-1000 mb) as does the Northern Hemisphere. The AZ(300-500 mb) 
results are presented for comparison with the meridional temperature 
gradients of Figure 4.1. 



WINTER 

SUMMER 

_ . 


January 

1970 

J anuary 
1970 

July 

1969 

July 

1969 

July 

1969 

January 

1970 


NMC 

N.H. 

SIRS 

N.H. 

SIRS 

S.H. 

NMC 

N.H. 

SIRS 

N.H. 

SIRS 

S.H. 


41.2 

38.4 

43.7 

16.3 

D 

34.9 

■ 


10.0 

16.0 

4.5 


9.1 


11.1 

4.0 

1 

3.2 

6.6 

H 

2.0 


2.0 


1.0 

1.8 

Hi 

0.6 


Table 4.1. Zonal and eddy available potential energy in the 

Northern and Southern Hemispheres for a winter and a 
summer month. All values are in units of 10 s 
joules m“ 2 . 




56 


Newell et al. (1970) present calculations of components pf thf 
atmosphere’s energy cycle in both hemispheres. The calculations of the 
zonal components are based on climatological data. The results of the 
present study in terms of AZ (100-1000 mb) confirm, at least qualita- 
tively, the hemispheric differences noted in the Newell study. The 
results are difficult to compare quantitatively because of differences 
in the computational technique. The calculations by Newell et al. cover 
the area 0°-90° latitude compared to the 20° -80° latitude span used in 
this study. In general, Newell's values for AZ(100-1000 mb) are larger. 

Evaluation of the time-mean zonal available potential energy can 
be made quite accurately from monthly mean, or even climatological 
temperature fields. For example, calculations of AZ(100-1000 mb) for 
July 1969 and January 1970 made on a daily basis from the NMC data 
fields, then averaged over the month, are only one percent larger than 
the same quantity calculated from the monthly-mean data fields. Quite 
the opposite is true for eddy terms, such as the eddy available potential 
energy. Estimations of AE made from monthly-mean temperature fields 
represent only the contribution of standing eddies. To obtain accurate 
estimates of AE, the computations must be made on a daily basis. 

The last two lines of Table 4.1 contain the results of the AE 
evaluations. Unlike the AZ calculations, those for AE indicate a large 
underestimation of the eddy term by the SIRS-based calculations. The 
reasons for this underestimation are discussed in Chapter 3 (Section 
3.5). Because of the large underestimation, inter-hemispheric compari- 
sons are difficult. Although the results for the 100-1000 mb layer are 
also presented, this author considers the computations for the 
300-500 mb layer more reliable for comparisons. The 300-500 mb layer 
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is the highest that is still basically part of the troposphere. 
Therefore, it should still exhibit the characteristics of the tropo- 
sphere and yet be relatively free (compared to lower layers) of cloud- 
contamination problems. Very tentative conclusions based on 
AE (300-500 mb) are that, at least in the upper troposphere, in winter 
the Southern Hemisphere has greater AE than the Northern Hemisphere, 
and in summer the two hemispheres have approximately equal quantities 
of AE. 

In winter AE (100-1000 mb) is calculated to be larger in the 
Northern Hemisphere, while AE (300-500 mb) is larger in the Southern 
Hemisphere. This variation with height in the hemispheric difference 
is probably true. At low levels the greater land-sea contrast in the 
Northern Hemisphere produces greater east-west gradients of temperature 
than over the relatively homogeneous Southern Hemisphere. Based on 
computations in this study, in winter AE(700-1000 mb) in the Northern 
Hemisphere is more than twice as large as that for the Southern Hemi- 
sphere. 

The computed KZ and KE for both hemispheres and months are given 
in Table 4.2. The values for KZ (75-1000 mb) and KE (75-1000 mb) are 
computed from height fields from 100 mb to 700 mb. KZ (150-250 mb) and 
KE(lS0-250 mb) are based on just the 200 mb heights. In both seasons 
the KZ values are larger in the Southern Hemisphere, especially in the 
summer. The KE values in Table 4.2 show little, if any, hemispheric 
difference, at least for the values for 200 mb (last line in table). 

At 200 mb, although the KZ is greater in the Southern Hemisphere in both 
seasons, the KE is approximately the same. This indicates, tentatively, 
that the ratio of KZ to KE in the Southern Hemisphere is greater than 
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that in the Northern Hemisphere. This is probably true for both seasons, 
but especially true in summer. 


• 

WINTER 

SUMMER 


January 

1970 

January 

1970 

July 

1969 

July 

1969 

July 

1969 

J anuary 
1970 




SIRS 

S.H. 

NMC 

N.H, 

SIRS 

N.H. 

SIRS 

S.H. 

KZ (75-1000 mb) 

18.2 

14.2 

17.4 

2.9 

2.5 

10.3 

KZ(150-2S0 mb) 

1 

5.3 

4.2 

5.2 

0.9 

0.8 

2.6 

KE (75-1000 mb) 

10.0 

6.0 

4.9 

5.1 

! 3.6 

3.5 

KE (150-250 mb) 

1.8 

i 

0.9 

0.9 j 

1.2 

0.6 

0.7 


Table 4.2. Zonal and eddy kinetic energy in the Northern and 

Southern Hemispheres for a winteT and a summer month. 
All values are in units of 10 5 joules m -2 . 


4.1.4 Standing and transient eddy contributions 

The eddy available potential energy (AE) can be thought of as 

consisting of the energy in standing eddies plus the energy in transient 

eddies. The percentage of AE contained in standing or stationary eddies 

is given in Table 4.3. The energy of the standing eddies is computed 

from the mean temperature structure for the month. Unlike the values 

for AE in Table 4.1, the percentages in the SIRS and NMC columns in 

* 

Table 4.3 for the Northern Hemisphere are very much alike. Thus, at 
least with regard to AE, it appears that the underestimation by the 
SIRS-based calculations is essentially the same with regard to both 
standing and transient waves. 

The percentages for the Southern Hemisphere in Table 4.3 are 
substantially lower than those for the Northern Hemisphere. The 
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inter-hemispheric difference is largest in summer with the Southern 
Hemisphere percentages being about one half the Northern Hemisphere 
values. In winter in the Northern Hemisphere, there is a significant 
contrast between the percentages for AE (100-1000 mb) and those for 
AE (300-500 mb). The higher percentages for the 100-1000 mb layer occur 
because of the contribution from the lowest layer (700-1000 mb). The 
land-sea contrast produces a large percentage (60%) of standing to total 
AE in that layer. No such difference in the percentages between 
AE (100-1000 mb) and AE (300-500 mb) is obvious in the Southern Hemisphere 
in either season, or in the Northern Hemisphere in summer. 



WINTER 


SUMMER 


January 

1970 

January 

1970 

July 

1969 


July 

1969 

July 

1969 

January 

1970 


NMC 

N.H. 

SIRS 

N.H. 

SIRS 

S.H. 


NMC 

N.H. 

SIRS 

N.H. 

SIRS 

S.H. 

2 

52% 

55% 

34% 


58% 

55% 

25% 


| 40% 

i 

42% 

31% 


54% 

58% 

23% 


Table 4.3. Percentage of eddy available potential energy (AE) 
contained in standing eddies. 


The percentages in Table 4.3 seem to indicate that standing eddies 
are more dominant relative to transient eddies in the Northern Hemisphere 
than in the Southern Hemisphere. This conclusion agrees with observa- 
tions by other authors (see Section 2.5). However, since the present 
study is based on only one month of data for each season, the reliability 
of this result will be examined closer. The analysis will focus on the 
winter comparison of AE (300-500 mb), where the smallest inter-hemispheric 
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difference occurs. It will be shown that there is a high probability 
that the difference in statistics is due to an actual difference between 
the hemispheres rather than to sampling limitations. 

To evaluate the probability outlined above, many years of calcula- 
tions of the percentages in Table 4.3 would be welcome. Such statistics 
do not exist. However, Ellis (1974) has calculated various general 
circulation parameters including AE (200-850 mb) for nine years 
(1963-1971). These statistics are examined in order to estimate the 
variances of the percentages in Table 4.3. For the AE (200-850 mb) 
parameter for nine Januarys, the values of Ellis show a sample mean of 
8.78 x 10 5 joules m" 2 and a standard deviation (s) of 0.79 x 10 5 
joules m" 2 . The ratio of standard deviation s to the mean is .09. It 
is assumed that this ratio of .09 holds for AE (300-500 mb) as wellj 
therefore, the standard deviation of this parameter is estimated to be 
0.06 x 10 5 and 0.09 x 10 5 joules m~ 2 for the Northern and Southern 
Hemispheres, respectively. In this exercise values of AE (300-500 mb) to 
two decimal places are used (the values in Table 4.1 are rounded-off 
to one place) . 

The term of interest here, however, is the ratio of AE (300-500 mb) 
due to standing eddies to the total AE (300-500 mb). This ratio multi- 
plied by 100 is identified by R. That is, 

= AE (standing eddies! x iqq, 4.9 

AE (total) ’ 


_ AE(standin 

or ~ AE (standing eddies) + 


eddies) 

(transient eddies) 


x 100. 


4.10 
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In the following paragraph let X represent the AE (standing eddies) and 
Y represent AE (total). Therefore, 

R = - x 100. 4.11 

Differentiating equation 4.11 gives 


dR 



x 100. 


4.12 


To estimate the standard deviation of R, o , the differentials in 
equation 4.12 are replaced by finite differences corresponding to the 
standard deviations of the variables. 


o 


R 


5 



x 100 


4.13 


Here X and Y represent typical values or means. It is assumed that 
deviations in X and Y are always of the same sign. All the variables 
on the right side of equation 4,13 are known except a . Equation 4.13 

X 

will be solved with a = 0 and with 

x 


a a 

y = 0,5 -y . 4.14 

Using the values in Table 4.4, is calculated in both hemispheres to 
be 3 is = 0 and 2 if is given by equation 4.14. Henceforth, the 
value of * 2.5 is used. 

It is assumed that the distribution being sampled from is Normal 
with an unknown mean y and standard deviation = 2.5. The highest 
probability that the respective Northern and Southern Hemisphere values 
of R (42 and 31) came from the same distribution occurs if y is halfway 
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between the two values 42 and 31; i.e., y = 36.5, The two values, 42 
and 31 are then both more than 2o away from y = 36.5. Accordingly, the 
probability that the two hemispheres are actually alike despite the 
differences in the percentages in Table 4,3 is low. Conversely, the 
probability is high that there is an actual disparity between the 
hemispheres . 



N. H, 

S. H. 

X 

.26 x 10 5 

.31 x 10 5 

Y 

.67 x 10 5 

1.01 x 10 5 

a 

.06 X 10 5 

.09 x 10 5 

y 



o x (by eq. 4.14) 

.01 x 10 5 

.01 x 10 5 

0r(o x - 0) 

3 

3 

j OrC",, by eq- 4.i4) 

2 

2 


Table 4.4. Values used in estimation of o R 


Table 4.5 shows the percentage of KB attributable to standing 
eddies. These values indicate the same hemispheric differences as are 
observed in the values for AE in Table 4.2. The percentages based on 
the 200 mb level Clast line of Table 4.5) show that the Southern Hemi- 
sphere percentages at tropopause level are less than those in the 
Northern Hemisphere. The Northern Hemisphere values are about the same 
in the two months, while in the Southern Hemisphere the summer values 
are lower. The same sort of trends are noted in the percentages for 
AE given in Table 4.3, The reason for this difference in seasonal 
variation of standing eddy importance in the two hemispheres lies with 
the relative strength of the summer monsoon in each hemisphere. The 
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presence of a very strong summer monsoon circulation in the Northern 
Hemisphere maintains, or even increases, the percentages shown in 
Tables 4.3 and 4.5. In the Southern Hemisphere, the summer monsoon 
circulation is much weaker, and the area-averaged percentages in the 
tables decrease from winter to summer. 



WINTER 

SUMMER 


J anuary 
1970 

J anuary 
1970 

July 

1969 

July 

1969 

July 

1969 

January 

1970 


NMC 

N.H. 

SIRS 

N.H. 

SIRS 

S.H. 

mjn 

SIRS 

N.H. 

SIRS 

S.H. 

KE(7S-1000 mb) 

1 31% 

33% 

23% 

28% 

34% 

17% 

KE (150-250 mb) 

32% 

31% 

24% 

30% 

31% 

18% 


Table 4.5. Percentage of eddy kinetic energy contained in 
standing eddies. 


4.2 Temperature variance and eddy kinetic energy as a function of 
latitude 

4.2.1 Temperature variance 

The contribution to eddy available potential energy from a specific 
layer, e.g. AE (300-500 mb), is a function of the variance of the tempera- 
ture around latitude circles weighted by cosif>. This variance as a 
function of latitude in the 300-500 mb layer is shown in Figure 4.4. In 
the Northern Hemisphere, curves are shown for both the NMC-based and 
SIRS-based calculations. There is a definite lack of agreement between 
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N.H. NMC 

Jan. 1970 

N.H. SIRS 

Jan. 1970 

— S.H. SIRS 
July 1969 


* 



N.H. NMC 

July 1969 

N.H. SIRS 

July 1969 

— S.H. SIRS 

Jan. 1970 


<t> 


Fig. 4.4 Variance of temperature around latitude circles for the 
layer 300-500 mb in both hemispheres; (a) winter, 

(b) summer. 
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the NMC and SIRS curves beyond the typical underestimation. In winter 
in the Northern Hemisphere, the SIRS curve indicates the maximum 
variance at 55^, while the NMC data indicates the peak at 45°N. The 
Southern Hemisphere temperature variance pattern in winter is quite 
different from that in the Northern Hemisphere. The variance appears 
to peak at a lower latitude in the Southern Hemisphere, but the reli- 
ability of this conclusion is questionable in view of the difference 
between the NMC and SIRS curves for the Northern Hemisphere. However, 
the Southern Hemisphere curve shows a secondary maximum in high latitudes 
not evident at all on the Northern Hemisphere curves. This high-latitude 
maximum is no doubt related to presence of the Antarctic continent. The 
relative contribution of standing and transient eddies in this maximum 
will be discussed in a following paragraph. 

The curves for summer (Figure 4.4b) show that the SIRS distribution 
does not adequately describe the features shown in the NMC curve for the 
Northern Hemisphere. The NMC curve shows two maxima, one at about 30°N 
and another at 55°N while the SIRS curve indicates the low latitude 
maximum, but not the peak at 55°N, The Southern Hemisphere curve shows 
only one maximum, at 55°S. The low latitude maximum in the Northern 
Hemisphere is related to the summer monsoon over Asia, and the maximum 
at 55°N is associated with migratory and stationary waves in the main 
belt of strong westerlies which has moved poleward from winter to 
summer. In the Southern Hemisphere in summer there is no temperature 
variance peak in low latitudes according to the SIRS-derived data. 

The 300-500 mb layer has been singled out in this section because 
it is the highest layer of the eight thickness layers obtained from the 
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satellite data that still reflects tropospheric temperature character- 
istics. The agreement between NMC and SIRS-based calculations increases 
with height. For example. Figure 4.5 shows the mean temperature vari- 
ance in the 200-300 mb layer. Agreement between NMC and SIRS curves in 
the Northern Hemisphere is substantially better than in the 300-500 mb 
layer just below. However, the variances in this upper layer, especially 
in high latitudes, are beginning to reflect occurrences in the strato- 
sphere instead of in the troposphere. 

The ratio of standing eddy temperature variance to total variance 
as a function of latitude indicates the latitudinal variability of the 
relative importance of standing and transient eddies. The distribution 
of this parameter for the 300-500 mb layer is given in Figure 4.6. 
Although there is rather poor agreement between the NMC and SIRS temper- 
ature variance curves in Figure 4.4, there is excellent agreement with 
regard to the ratio in Figure 4.6. The percentages in the second row 
of Table 4.3 in Section 4.1.4 are essentially area or cosine-weighted 
averages of the ratio in Figure 4.6. In both seasons there are signif- 
icant differences between the hemispheres. In winter both hemispheres 
show a relatively high value at 20° latitude, where the total variance 
is small, and a relative minimum at 25° or 30° latitude. Poleward of 
that latitude, the patterns are quite different. North of 30^, there 
is a single maximum in mid-latitudes (50°-55°N), while the Southern 
Hemisphere shows a double maximum with a peak at 40° S and a general area 
of high values at, and south of, 65°S, The areas of maximum relative 
contribution by standing eddies coincide generally with the maxima of 
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N.H. NMC 
Jon. 1970 

N.H. SIRS 
Jan. 1970 

S.H. SIRS 
July 1969 





N.H. NMC 
July 1969 

N.H. SIRS 
July 1969 

S.H. SIRS 
Jan. 1970 


Fig. 4.S 


Variance of temperature around latitude circles for the 
layer 200-300 mb in both hemispheres; (a) winter, 

(b) summer. 
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N.H. NMC 
Jon. 1970 

N.H. SIRS 
Jan. 1970 

S.H. SIRS 
July 1969 
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N.H. NMC 
July 1969 

N.H SIRS 
July 1969 

S.H. SIRS 
Jan. 1970 


Fig. 4.6 Ratio of temperature variance due to standing eddies to 
total variance for the layer 300-500 mb; (a) winter, 
(b) summer. 
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total variance (Figure 4.4a). That is, large temperature variance is 
associated with a relatively large percentage contribution by standing 
eddies . 

In summer (Figure 4.6b) the latitudinal distributions in each 
hemisphere are similar with regard to the location of relative and 
absolute maxima and minima, but not with regard to magnitudes. At all 
latitudes the Northern Hemisphere values exceed those of the Southern 
Hemisphere. There has also been a large change in the patterns from the 
winter to the summer month. Both hemispheres display an absolute maxi- 
mum in the region of 25-30° with the Northern Hemisphere values exceeding 
those of the Southern Hemisphere. This maximum is related, at least in 
the Northern Hemisphere, to the summer monsoon. In the Southern Hemi- 
sphere, the total variance is small (Figure 4.4b) and the percentage 
contribution of standing eddies is also relatively small. These two 
facts merely indicate the relative weakness of the Southern Hemisphere 
summer monsoonal regime. Both hemispheres have a minimum at 45° and 
relatively high values poleward of about 55°. The Northern Hemisphere 
values are significantly higher than those in the Southern Hemisphere. 
These maxima around 60° in Figure 4.6b are associated with the maxima of 
total variance at 55° in Figure 4.4. 

4.2.2 Eddy kinetic energy at 200 mb 

The distributions of [u]. and, therefore, the zonal kinetic 

A,t 

energy are well known in the Southern Hemisphere because they can be 
calculated accurately from climatological data. However, the distribu- 
tion in space of eddy kinetic energy, especially at upper tropospheric 
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and tropopause levels, is not well known. Kao and Jenne and Sagendorf 
(1970) and van Loon (1972c) use the I GY 500 mb height fields to examine 
the latitudinal distributions of eddy kinetic energy in the Southern 
Hemisphere. Both studies show mid-latitude peaks (~ 40-50°S) in summer 
and winter, similar to what is found in the Northern Hemisphere (see 
Oort and Rasmusson, 1971) at 500 mb. The distribution of eddy kinetic 
energy at 200 mb is, however, quite different from that at 500 mb. 

Eddy kinetic energy at 200 mb as a function of latitude is shown 
for winter and summer in Figures 4.7 and 4.8, respectively. The solid 
curves represent the total eddy kinetic energy, and the dashed curves 
represent the standing eddy contribution. The difference between the 
two curves is the transient eddy kinetic energy. Parts a and b of each 
diagram show the Northern Hemisphere values calculated from NMC data and 
from the SIRS-based structure, respectively. The NMC curves for winter 
(Figure 4.7a) indicate one distinctive maximum at 32.5°N for both the 
total and standing eddy kinetic energy. This agrees well with Oort and 
Rasmusson (1971), whose tables show a maximum at 35°N. The SIRS curves 
(Figure 4.7b) also show the maxima at 32.5°N, but, in addition, show 
secondary maxima at 57.5°N. The appearance of this second peak is 
totally related to the standing eddy maximum. This erroneous peak in 
the standing eddy contribution at 57.5°N is probably related to the 
problems in the distribution of satellite data which produced the errors 
in the [zK curve at 200 mb (Figure 4.2a) and in the [T] curve for 

A j t 

the 300-500 mb layer (Figure 4.1a) in January in the Northern Hemisphere 
at the same latitude. This data distribution problem does not occur in 
the other month or the other hemisphere. 
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Fig. 4.7 Distribution with latitude of eddy kinetic energy at 200 mb 
in both hemispheres during winter. Total eddy kinetic 
energy (solid lines), standing eddy kinetic energy (dashed 
lines). 
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Fig. 4.8 Distribution with latitude of eddy kinetic energy at 200 mb 
in both hemispheres during summer. Total eddy kinetic 
energy (solid lines), standing eddy kinetic energy (dashed 
lines) . 
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The Southern Hemisphere winter (July 1969) curves are shown in 
Figure 4.7c. The peaks in the total eddy kinetic energy and in the 
standing eddy contribution occur at 27.5°S, at a slightly lower latitude 
than those in the Northern Hemisphere. This observed difference may 
very well just be a function of the small sample. The magnitudes of 
the peaks in the two hemispheres in winter, as calculated from the 
satellite-based structure, are nearly the same. The Southern Hemisphere 
curve also shows a secondary maximum in high latitudes, poleward of 70°S 
(Figure 4.7c). No similar peak exists in the Northern Hemisphere. This 
peak in the eddy kinetic energy is composed of maxima in both the 
standing and transient eddy components and also shows up in calculations 
for lower layers, including 500 mb. However, the calculations of Kao, 
Jenne and Sagendorf (1970) do not show this peak in the IGY data at 
500 mb, although their computations do not extend farther poleward than 
70°S. The magnitude of the high latitude peak in Figure 4.7c may also 
be larger than what it normally is. The evidence for this are the large 
values of [uL from 60° to 75°S (Figure 4.3a), relative to the clima- 

A,t 

tological values at that latitude. If the eddy kinetic energy is 
associated with this jet stream, then the values in Figure 4.7c may be 
an overestimate of normal conditions. Therefore, the existence under 
normal conditions of an eddy kinetic energy maximum at high latitudes 
is undecided. 

At the latitude of peak eddy kinetic energy in both hemispheres in 
winter, the standing eddies contribute significantly. In the Northern 
Hemisphere at 32.5°N, the standing eddies percentages are 39% and 41% 
for the NMC and SIRS-based calculations, respectively. At 27.5°S the 
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percentage as determined from the SIRS-based calculations is 30%. These 
values are larger than the percentages averaged over the hemisphere, 
presented in Table 4.5. This means that there is a positive correlation 
between the eddy kinetic energy and the percentage contribution of 
standing eddies. At 500 mb, Kao, Jenne and Sagendorf (1970) show the 
peak eddy kinetic energy to be at about 50°S in winter. The percentage 
contribution of standing eddies there is 12% according to their calcula- 
tions. This percentage is much lower than the 30% at 200 mb calculated 
here. This disparity may be related to the difference in level used, a 
difference in the latitude location of the peak, or to an underestimation 
of the standing eddy contribution by Kao et al. However, at 200 mb and 
50°S the percentage is still about 30%, using the statistics of the 
current study. Another reason for the variation could be that, as far 
as can be determined, Kao's standing eddy contribution is computed from 
fields averaged over three months rather than the one month used in this 
study. This reduces the percentage contribution of the standing eddies. 
However, it seems unlikely that this difference in averaging time could 
account for all the difference between 30% and 12%. It therefore seems 
possible that studies based on the IGY data fields have underestimated 
the role of standing eddies. 

The distribution of 200 mb eddy kinetic energy in summer is shown 
in Figure 4.8. The peak has shifted poleward to about 45° in both 
hemispheres with the Southern Hemisphere magnitude being substantially 
larger than the Northern Hemisphere one, as calculated from SIRS-based 
structure in both hemispheres (Figures 4.8b and c) . 



4.3 Eddy heat transport 

4.3.1 Distribution of eddy heat transport 

The computations of eddy heat flux in the 300-500 mb layer reveal 
significant differences between the hemispheres. The fluxes are calcu- 
lated using the 300-500 mb thickness information and the 300 mb and 
500 mb height fields. Geostrophic wind speeds are used throughout. The 
300-S00 mb layer represents the upper troposphere. Although the relative 
contribution of this layer to the vertically integrated eddy heat flux 
is small, the layer should possess characteristics of the tropospheric 
eddy flux fields. 

Figure 4.9 shows the results of the eddy heat flux computations for 
the winter month in each hemisphere. The solid lines represent the 
total eddy flux, [(v)^(T)^]^ t , and the dashed lines represent just the 
standing eddy flux, [ ( [v] ) ^ ( [T] . Parts b and c of the figure have 

different ordinate scales than part a. Figures 4.9a and 4.9b are the 
result of the Northern Hemisphere computations based on NMC and SIRS- 
derived data, respectively. Both diagrams show a maximum in mid- 
latitudes, although not at exactly the same latitude. The underestima- 
tion of the total eddy flux is about 50% at the latitudes of the maxima. 
The curves for the standing eddies also show a miximum in mid-latitudes. 

In addition, it is evident from Figures 4,9a and b that the 
transient eddy flux is underestimated more than the standing eddy flux. 
This is quite different from the situation in the second section of this 
chapter where the relative importance of standing and transient eddies 
in the variance of temperature around latitude circles was well por- 
trayed by the SIRS-based calculations (see Figure 4.6). The reason for 
this difference lies in the size of the eddies involved. The variance 
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Fig. 4.9 Wintertime, total poleward eddy flux of heat, ± [ (v) A (T) A h 

(solid lines), and poleward eddy flux of heat due to standing 
eddies, ± [ ( [v] t ) A ( [T] t ) A ] A (dashed lines) for the layer 
300-500 mb in both hemispheres. 
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of temperature around latitude circles in the upper troposphere is 
dominated by large eddies of hemispheric wave numbers 1 through 3 
(Saltzman and Fleisher, 1960) . The eddy flux of heat in mid-latitudes 
is dependent for the most part on large, stationary waves (wave number 
~ 3) and shorter (wave number ~ 6), transient waves (Kao and Sagendorf, 
1970). The shorter, transient waves are not picked up as well by the 
SIRS-based structure as are the larger, stationary or near- stationary 
waves. Therefore, the flux due to stationary waves is better estimated 
by the SIRS-based calculations than that due to transient waves. 

There are striking variations between the curves in Figure 4.9 for 
the Northern and Southern Hemispheres. First, the maximum eddy flux in 
the Southern Hemisphere (Figure 4.9c) occurs at 3S 0 S, farther equatOTward 
than either calculation in the Northern Hemisphere (Figures 4.9a and b) . 
However, because the SIRS-based calculation in the Northern Hemisphere 
(Figure 4.9b) does not correctly identify the latitude of maximum flux, 
any conclusion as to the relative positions of the latitude of maximum 
flux in the two hemispheres would be premature. This is especially true 
since the calculation is based on only one month of data. Other calcu- 
lations which show the seasonal variation of eddy heat flux (Newell, 
et al . , 1970; Kao, Jenne and Sagendorf, 1971) in the Southern Hemisphere 
indicate a peak at about 50 6 S in winter. 

The standing eddies show two peaks of poleward heat flux in the 
Southern Hemisphere, one at 35 o -40°S and another at 70°S. In mid- 
latitudes the flux curve has a relative minimum. This distribution 
generally parallels that of the percentage of temperature variance due 
to standing eddies given in Figure 4* 6a. Poleward of 65°S the heat flux 
is dominated by the standing eddies. In mid-latitudes the percentage 
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of the total eddy flux attributable to standing eddies is much less in 
the Southern Hemisphere than in the Northern Hemisphere. At the 
latitude of peak eddy flux in each diagram, the percentage due to 
standing eddies is shown in Table 4.6, 



N.H. 

S.H. 

January 1970 
NMC 

January 1970 
SIRS 

July 1969 
SIRS 

Standing Eddy 
Percentage 

54% 

67% 

24% 


Table 4.6. Percentage of eddy flux of heat in the 300-500 mb 

layer attributable to standing eddies at latitude of 
peak eddy flux. 


The Southern Hemisphere percentage in Table 4.6 is probably somewhat 
overestimated judging by the Northern Hemisphere NMC and SIRS percent- 
ages. Kao, Jenne and Sagendorf (1971) present heat flux statistics for 
winter, based on analyzed 500 mb height and 500-1000 mb thickness fields 
from the IGY, which give a percentage of 12% at 30°S. This value is 
much less than the one calculated in this study. The difference cannot 
be totally attributed to the overestimation by the SIRS calculation. 

The resolution of these differences pertaining to the relative importance 
of standing eddies in the 30°-40°S latitude band must await further 
study. Kao, Jenne and Sagendorf (1971) do indicate the standing eddy 
peak at about 30°S with a sharp decrease to 50°S, but do not indicate 
the high latitude maximum that is evident in Figure 4.9c, 

In summer (Figure 4.10) the standing eddies in each hemisphere are 
less important in mid- latitudes than in winter. The total eddy flux is 
less in summer in both hemispheres, but the seasonal change is much 
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Fig. 4.10 Summertime, total poleward eddy flux of heat, ± [(v)^(T), ^ 
(solid lines), and poleward eddy flux of heat due to standiftg 
eddies, ± [ ( [v] t ) x ( [T] ) x ] x (dashed lines) for the layer 
300-500 mb in both hemispheres. 
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larger in the Northern Hemisphere. The latitude of peak eddy flux has 
moved poleward to about 60° in each hemisphere. The Northern Hemisphere 
SIRS curve (Figure 4.10b) does not show the sharp drop-off with increas- 
ing latitude poleward of 60°N that is evident in the NMC curve (Figure 
4.1Ca). This discrepancy is tied in part to the failure of the SIRS- 
based calculations to correctly identify the large negative fluxes of 
the standing eddies north of 70°N. 

4.3.2 Mid-latitude magnitude 

Robinson (1970) has concluded, on the basis of calculations using 
primarily IGY station data, that the vertically integrated total eddy 
flux of heat in mid-latitudes is much greater (by a factor of two) on 
an annual average in the Northern Hemisphere than in the Southern 
Hemisphere. At 50° latitude at 400 mb, a similar-sized hemispheric 
difference is obtained from the tables and charts of Robinson. Newton 
(1972) concludes that Robinson's Southern Hemisphere mid-latitude 
fluxes are much too low. Instead he states that the magnitude of the 
eddy flux of heat at latitude 50° (the approximate latitude of maximum 
eddy flux) is about the same in both hemispheres. This conclusion is 
reached as a result of energy budget calculations. The eddy heat flux 
calculations in the present study offer an opportunity to compare the 
two hemispheres and possibly determine the relative strengths of the 
Northern and Southern Hemisphere's eddy fluxes. 

The inter-hemispheric comparison is made at latitude 50°, the 
approximate latitude of the annual maximum. Because of the contrast 
between the hemispheres in relative strength of standing eddies and the 
difference in the ability of S IRS-derived structure to determine 
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correctly the standing and transient eddies, the values in parts b and 
c of Figures 4.9 and 4.10 cannot be immediately used to establish the 
comparison. Therefore, the following procedure is adopted. The total 
eddy flux is divided in each hemisphere and each month into the standing 
and transient eddy portions . In the Northern Hemisphere the NMC and 
SIRS-based values are compared. The ratio of NMC to SIRS-based values 
gives a "correction factor" which when multiplied by the Northern 
Hemisphere SIRS values results, of course, in the NMC or assumed correct 
fluxes. The major assumption in this exercise is that the Northern 
Hemisphere "correction factors" apply in the Southern Hemisphere. 
Separate factors are determined for standing and transient eddies. For 
January 1970 in the Northern Hemisphere, the transient and standing eddy 
"correction factors" are 3.6 and 1.2, respectively. In July the tran- 
sient eddy factor is 4.1. There is no standing eddy correction term for 
July because the standing eddy flux at 50°N is negligible. The two 
transient eddy "correction factors" are averaged, and the mean value of 
3.8 is used in the Southern Hemisphere for both months. Applying the 
derived "correction factors" to the Southern Hemisphere fluxes, adding 
the standing flux to the transient flux, and averaging over January and 
July results in the value in the third column of Table 4.7. The 
adjusted, or corrected, value for the Southern Hemisphere is smaller, 
but only 14% smaller than the NMC, Northern Hemisphere value. The 
Southern Hemisphere corrected SIRS value is also approximately twice as 
large as the annual average value at 400 mb from Robinson (1970). The 
first entry in the table is from Oort and Rasmusson (1971) at 400 mb. 

It is smaller than the other two values. This may be because the NMC- 
based calculations overestimate the flux, or that 1969-1970 was a 
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period of above average flux. Table 4.7 does indicate, however, that the 
Southern Hemisphere flux is closer to being equal to the Northern 
Hemisphere value than to being one half of the Northern Hemisphere value. 
This conclusion agrees with Newton (1972) and disagrees with Robinson 
(1970). 



- . ... - - 

N.H. 

Climate 

N.H. 

NMC 

S.H. 

SIRS 

S.H. 

ROBINSON 

Upper Troposphere 
Total eddy heat flux 

7.0 

9.9 

8.5 

4.1 


Table 4.7. Upper troposphere "corrected" total eddy heat flux 

[(v)-^(T)^]^ t at 50° latitude averaged over January and 
July. N.H. ’ climate value from Oort and Rasmusson (1971) 
at 400 mb. Robinson (1970) S. H. value an annual mean 
at 400 mb. Units: °K m sec -1 . 

4.3.3 Longitudinal distributions 

The last few paragraphs have dealt mainly with the latitudinal 
distribution of eddy heat flux. The longitudinal variation at a partic- 
ular latitude is also important. Haines and Winston (1963) discuss the 
longitudinal distribution of eddy heat flux in the Northern Hemisphere 
for the layer 500-850 mb. At 45°N during the winter season (October- 
March) , there are three bands of longitude which contribute greatly to 
the longitudinally-averaged poleward flux of sensible heat. The three 
zones, obvious in Figure 4.11, are associated with the location of mean 
ridge and trough positions (Haines and Winston, 1963). The narrow band 
located at about 120°E is on the west side of the mean trough off the 
Asian coast. The peak in the eastern Pacific Ocean is situated on the 
east side of this trough. The general maximum between 20°W and 100°W 
is associated with the trough located at the east coast of North 




Fig. 4.11 Averages for four consecutive cold seasons (October-March) 
of contributions in each 5° longitude zone to poleward eddy 
heat transport across 45°N (from Haines and Winston, 1963). 
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America, with the peak values at about 40-50°W occurring on the east 
side of the trough. The relation of the peak locations in Figure 4,11 
to the standing eddy pattern in the troposphere does not mean, however, 
that the fluxes in the figure are totally due to the standing eddies. 

At 45 °N in winter, the eddy heat flux is divided about equally between 
transient and standing eddies. However, the transient eddies have 
preferred longitudes of greatest amplitude. 

Because the Southern Hemisphere standing eddies are relatively less 
important, it would be interesting to determine if there are still 
preferred longitudes of eddy heat flux in the Southern Hemisphere. For 
this reason calculations of the eddy flux of heat in the 300-500 mb 
layer for every 10° of longitude were carried out at 45°N and 45°S, at 
45°N using both NMC and SIRS-based structure. Figures 4.12a and b show 
the results for the Northern Hemisphere. In January the NMC curve 
clearly shows the three areas of poleward heat flux also observed by 
Haines and Winston (1963) for the lower troposphere. The SIRS curve for 
January (Figure 4.12a) indicates the strong peak at about 120°E, but 
only slightly indicates the other two maxima. In summer (Figure 4.12b) 
the patterns are much weaker (the ordinate scales are different in 
Figures 4.12a and b) and less cohesive. However, the SIRS curves do 
indicate the longitudes of large poleward fluxes. 

The Northern Hemisphere NMC total flux curve for January 1970 is 
reproduced in Figure 4.13 as a solid line. Also shown is the flux 
associated with the standing eddies (calculated from January 1970 mean 
temperature and height fields). Averaged around the latitude circle, 
the standing eddies contribute about 50% to the total eddy flux at 45°N 
(see Figure 4.6a). They obviously play an important role in the 
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Longitude (°E) 


Fig. 4.13 Longitudinal distribution of total eddy heat flux and 

standing eddy heat flux for the layer 300-500 mb at 45°N 
in January 1970, 
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longitudinal variations. The peak eddy flux at 125°E is due almost 
entirely to standing eddies. The other two regions of large poleward 
flux show a very pronounced east-to-west variation in the importance of 
standing eddies. Both regions (one from 200 o -240°E, the other between 
270° and 360°E) show the standing eddies dominating in the eastern part 
of the region, and transient waves being relatively more important in 
the western portion of each region. The contribution by transient eddies 
is given by the difference between the solid and dashed curves. The 
same east-to-west variation is present for the region of large flux 
centered at 120°E J but is not as obvious from the diagram. This east- 
to-west variation of standing and transient eddy importance reflects the 
presence of travelling systems amplifying and slowing as they move east- 
ward until they reach maximum amplitude and become stationary. 

The longitudinal variations at 45°S calculated from the SIRS-based 
structure are shown in Figure 4.14 for July 1969 and January 1970. The 
45°S latitude circle is land-free, except for where it crosses the 
Southern tip of South America and New Zealand. In the winter month, 

July 1969, 22% of the total eddy flux is calculated to have been 
attributable to stationary eddies. This percentage is probably a slight 
overestimate, because tmasient eddies are underestimated more than 
standing eddies (see Section 4.3.1). The calculated percentage for 
January 1970 is 6%. Both months show longitudinal variations, the winter 
month showing more organization. In the winter month there are three 
peaks in the poleward eddy flux: one at 35°E, east of the axis of the 

African continent; a second at 17S°E, just east of New Zealand; and the 
third from 285°E to 325°E, across and to the east of South America. 

Peaks in the stationary eddy flux (dashed line in Figure 4.14a) occur 



(°K m sec -1 ) (°K m sec"‘) 

45° S, 300-500 mb 45° S, 300-500 mb 

Poleward Eddy Heat Flux Poleward Eddy Heat Flux 
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Standing Eddy Flux 
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in conjunction with two of the total eddy flux peaks. The South 
American peak is essentially comprised of just transient eddy flux. The 
African total eddy flux peak consists of about one half stationary and 
one half transient eddy flux. 

An interesting feature in Figure 4.14a is the displacement of the 
stationary eddy flux maximum to the east of the total eddy flux peak 
over New Zealand. The displacement is 20° longitude. Between 185° and 
235°E, the stationary eddy flux is greater than the total; therefore, 
the transient eddy flux is negative. 

In summer (Figure 4.14b) the standing eddies are nearly negligible, 
except at a few longitudes, and the pattern of total eddy flux is more 
confused. The very narrow peak at 245°E is the result of a very short- 
lived, but very strong transient system. The other peak areas are in 
the general locations of their winter counterparts. In the Southern 
Hemisphere, the peak values in summer are about two-thirds of the winter 
maximum values. In the Northern Hemisphere, the summer value is only 
about 18% of the winter value. This again points out the large differ- 
ence between the hemispheres in seasonal variation. 

4.4 Chapter summary 

The most significant results in this chapter are those pertaining 
to the eddy terms. The mean meridional gradients and parameters, such 
as AZ and KZ, based on them merely confirm previous results based on 
climatology (Sections 4.1.2 and 4.1.3). It should also be mentioned 
that the mean meridional gradients calculated from the satellite data 
are constrained to be similar to the actual ones because separate 
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regression equations are used for each of three latitude bands in each 
hemisphere (see Section 3.3). 

Calculations presented in this chapter indicate that standing 
eddies are much less dominant in the Southern Hemisphere than they are 
in the Northern Hemisphere. This is true for the temperature structure 
as revealed in the computations of AE (Section 4.1.4), the eddy kinetic 
energy (Section 4.1.4) and the eddy heat flux (Section 4.3.1). The 
satellite-based values of the percentage contribution of standing waves 
to AE and KE are also shown to agree very closely to the NMC-based 
quantities. With regard to the eddy heat flux the satellite-based 
percentage is significantly larger than the corresponding NMC value. 

The variance of temperature around latitude circles is divided into 
standing and transient portions in Section 4.2.1. The ratio of standing 
eddy temperature variance to total variance as a function of latitude 
describes the latitudinal variation of the relative importance of 
standing eddies. These statistics show that in winter, in mid- latitudes, 
standing eddies are much more dominant in the Northern Hemisphere than 
they are in the Southern Hemisphere. This agrees with the hemispheric- 
averaged percentages in Table 4,3 in Section 4.1.4, However, in high 
latitudes (>65°), the reverse is true with the Southern Hemisphere 
standing eddies making a larger contribution to the total variance than 
their Northern Hemisphere counterparts. In summer the standing eddy 
contribution is larger in the Northern Hemisphere at all latitudes, the 
greatest difference occurring in low latitudes (<40° ). At these lati- 
tudes in the Northern Hemisphere, the standing eddy contribution is very 
large and is associated with the presence of the Asian monsoon. The 
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summer monsoon regime in the Southern Hemisphere is much weaker and the 
standing eddy contribution in low latitudes is much less. 

The magnitude of the mid-latitude eddy heat flux averaged over 
January and July is shown to be approximately equal in the two hemi- 
spheres (Section 4.3,2). This is significantly different from a previous 
observational study based on IGY data (Robinson, 1970) which concluded 
that the Southern Hemisphere eddy flux is only one half as large as 
that in northern mid-latitudes. 

The eddy heat flux in southern mid-latitudes is also shown to have 
significant longitudinal variations (Section 4.3.2). In winter, when 
the pattern is best defined, there are three peaks of poleward heat 
flux, apparently associated with the three land areas, Australia, South 
America, and Africa. 



CHAPTER 5 


A COMPARISON BETWEEN HEMISPHERES OF 


THE STRUCTURE OF 


STRATOSPHERIC POLAR VORTICES 


5.1 Relation of vortex structure to breakdown characteristics 

5.1.1 Background and theory 

In Section 2.6 it was noted that the Northern Hemisphere strato- 
spheric vortex usually breaks down in midwinter, while the Southern 
Hemisphere vortex does not. A possible explanation for this disparity 
in breakdown climatologies is a difference in the stability of the 
vortices in the two hemispheres. Another possibility related to the 
upward propagation of Rossby waves is discussed in Section 2.6. In 
this section the stability of the vortices in both hemispheres will be 
examined. It is shown that the Southern Hemisphere vortex appears to 
be stable in midwinter, but not in spring when it is observed to break 
down. 

Considering a fluid with both horizontal and vertical shears, 

Chamey and Stem (1962), Pedlosky (1964), and Mahlman (1966) derive 

essentially identical stability criteria. Chamey and Stem (1962) and 

Mahlman (1966) apply their criteria to the stratospheric polar vortex. 

As derived in these studies, the necessary condition for instability is 

3 0 

that the latitude profile of potential vorticity, P = (- (Sg+f ) > must 

have a relative maximum or minimum. In other words, [(- g^Htg+f)] 
must be zero at some point in the latitude profile. In the expression 
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for the potential vorticity, 0 is the potential temperature, p is 
pressure, is the relative vorticity on an isentropic surface and f 
is the Coriolis parameter. 

Charney and Stern (1962) show latitude profiles of potential vor- 
ticity for a case study just before a breakdown of the vortex in the 
Northern Hemisphere* A relative maximum in P is located just poleward 
of the jet axis; therefore, the instability criterion is met* The 

latitude profiles of the components of the potential vorticity, the 

30 

static stability (- ~) and absolute vorticity (C Q +f), indicate that the 
distribution of static stability is of primary importance in meeting the 
instability criterion (Mahlman, 1966). While the absolute vorticity 
shows a continuous increase with increasing latitude, the static stabil- 
ity decreases sharply from the jet axis poleward. If the static 
stability had remained constant with latitude, the instability criterion 
would not be met. Thus the distribution of static stability with 
latitude is the determining factor in the vortex stability criterion. 
This is fortunate for the present study, because calculations of ? +f 
are not available. The vortex structure is examined solely with respect 
to the static stability. 

i 

5.1.2 Some observations of vortex structure 

The latitudinal distribution of static stability is quite different 
in the stratospheric polar vortices of the Northern and Southern Hemi- 
spheres. In midwinter, poleward of the axis of the stratospheric jet 
stream, the stability decreases toward the pole in the Northern Hemi- 
sphere. Figure 5.1 shows the distribution with latitude of static 
stability in the layer 50-100 mb. The solid line indicates the latitude 
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Fig. 5.1 Static stability (- ■— .) in the layer 50-100 mb as a function 
of latitude in the ^ Northern Hemisphere for January. 
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Fig, 5.2 Static stability (- in the layers 50-100 mb and 30-50 mb 
as a function of latitude in the Southern Hemisphere for July. 
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profile obtained from January mean cross-sections presented by Craig 
(1965), which are based on four years of data. The decrease of stability 
with latitude is evident, the value at 90°N being 10% less than that at 
70°N. Because vortex breakdowns often occur during January, the January 
mean thermal structure is not the best example of pre-breakdown condi- 
tions. The dashed line in Figure 5.1 displays an example calculated by 
Mahlman (1966) of the stability profile before the onset of a major 

.1 

breakdown in 1958. A sharp decrease with latitude is again present, 
with the value of static stability at 90 o N being about 19% less than 
that at 70°N. 

The midwinter structure for the Southern Hemisphere is given in 
Figures 5.2 and 5.3. The dashed line in Figure 5.2 gives conditions in 
the 50-100 mb layer in July 1957 from station data presented by Taylor 
(1961). The stability at 90°S is calculated from July 1957 mean temper- 
atures for station Amundsen-Scott . The point plotted at 69°S is an 
average of conditions at Cape Hallett (72°S, 170°E) and Wilkes (66°S, 
111°E). The static stability is nearly constant with latitude with onl} 
a 4% decrease from 70°S to 90°S. The profile for July 1969 for the 
layer 30-50 mb is given by the solid line in Figure 5.2 and is based on 
temperature maps presented by Labitzke and van Loon (1972) . Again fror 
70°S to 90°S there is a decrease of 4% in stability. 

Figure 5.3 shows the latitude profile of hydrostatic stability ii 
the Southern Hemisphere based primarily on temperatures obtained from 
the SIRS data. The point at 90°S is based on station Amundsen-Scott. 

The curve in Figure 5.3 is much like the profile given by the solid line 
in Figure 5.2, with the exception of the increase in static stabilit; 
from 80°S to 90°S. The lapse rate at 90°S is based on only three da s 
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Fig. 5.3 


3 0 

Static stability (- t— ) in the layer 40-75 mb as a function 
of latitude in the P Southern Hemisphere for July 1969. 
Values (except at 90°S) determined from SIRS data. 
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SOUTH LATITUDE 

Fig. 5.4 Absolute vorticity as a function of latitude at 30 mb in 
the Southern Hemisphere for July 1969, 
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of data at the higher level, therefore it is probably not an accurate 
measure of the monthly mean. However, it appears that there is no large 
decrease in high latitudes as there is in the Northern Hemisphere. 

The latitude range from 70° to 90° is of the most importance for 
the following reason. Figure 5.4 gives the latitude profile of absolute 
vorticity at 30 mb for July 1969 as calculated from the distribution 
with latitude of zonal wind given by Labitzke and van Loon (1972) . While 
the vorticity on a pressure surface is not the same as the vorticity on 
an isentropic surface, the author believes that the profile in Figure 5.4 
is characteristic of conditions occurring on isentropic surfaces in the 
mid-stratosphere. The absolute vorticity decreases (becomes larger 
negative) with latitude, with the maximum rate of decrease occurring at 
60°S at the location of the jet axis. Poleward of about 70°S, the 
absolute vorticity profile flattens, becoming nearly constant with 
latitude. It is in this latitude range that a sharp decrease of static 
stability with increasing latitude will produce an extremum in potential 
vorticity, the criterion for vortex instability. This appears to be the 
case in the Northern Hemisphere where a sharp decrease in static stabil- 
ity with latitude exists in high latitudes (see Figure 5.1). The 
Southern Hemisphere, however, does not exhibit the large decrease from 
70°S to 90°S (see Figures 5.2 and 5.3). Thus, while the instability 
criterion is met during midwinter in the Northern Hemisphere (Mahlman, 
1966; Chamey and Stern, 1962), it will tend not to be met during 
midwinter in the Southern Hemisphere. 

Another indication of the difference in the early winter (before 
Northern Hemisphere breakdown) thermal structure of the Arctic and 
Antarctic vortices is given in a comparison of satellite-observed 



98 


radiances over the two polar regions, discussed by Labitzke and van Loon 
(1972). Although the Southern Hemisphere pole is much colder in mid- 
winter at 30 mb, radiances in SIRS channel 8 (665 cm -1 ) and in channel 
A of the Selective Chopper Radiometer (SCR) indicate that the upper 
stratosphere over the Antarctic is warmer than over the Arctic. There- 
fore, the static stability would be larger over the South Pole than over 
the North Pole. At about 50° latitude the 30 mb temperature and the 
radiances indicate that conditions are much alike in the two hemispheres. 
Thus, the gradient of static stability (- -^) is dissimilar in the two 
hemispheres. 

The variation in static stability distributions between the 
Northern and Southern Hemispheres in winter is probably related to one 
of, or a combination of, the following factors: a difference between 

the hemispheres in the vertical distribution of radiative cooling rates 
in the polar -night region; a difference in the patterns of eddy heat 
fluxes; or a difference in mean meridional circulations. The last two 
factors are investigated in Section 5.2. 

The distribution of total radiative heating in the global atmosphere 
is discussed by Dopplick (1972) on a season-by-season basis. His cal- 
culations reveal a large difference between the hemispheres in the 
magnitude of the total radiative cooling in the polar stratosphere during 
winter. The total cooling rate is substantially larger in the northern 
polar region. Dopplick relates this to the much lower Southern Hemi- 
sphere temperatures in this region of the lower stratosphere. The 
cooling rate also increases with height in the polar regions of both 
hemispheres; therefore, radiation acting alone will result in static 
stability decreasing with time. The distribution with latitude of this 
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destabilizing effect also varies between Northern and Southern Hemi- 
sphere high latitudes, At 70° and 80° latitude (Dopplick's computations 
extend to only 80°) the destabilizing effect of radiation is the same in 
the Northern Hemisphere, while in the Southern Hemisphere it decreases 
toward the higher latitude. Therefore, the Northern Hemisphere distri- 
bution will tend to produce a temperature structure with static stability 

/ 

constant with latitude, while static stability increasing toward the 
pole would result in southern latitudes. Thus the radiation distribution 
itself does not seem to produce the structure observed, although the 
inter-hemispheric difference in the distribution of radiative cooling in 
association with the other effects may be important in the total expla- 
nation of the observed variations in vortex structure. 


5.1.3 Winter- to-spring change in the Southern Hemisphere 
In spring as the sun returns to high latitudes in the Southern 
Hemisphere, changes in the stratospheric lapse rate occur in such a way 
as to produce a larger latitude gradient of static stability. The 

heating due to the return of the sun and the presence of ozone results 

3 0 

in increases in (- -t—) , with greater increases occurring at lower 

9p 

latitudes because of the greater insolation. Figure 5.5 shows the 
changes in stability in the 50-100 mb layer which take place between 
June and October, based on monthly mean data from Amundsen-Scott, Cape 
Hallett and Wilkes stations in 1957. The data for Cape Hallett and 
Wilkes is again averaged and plotted at latitude 69°S, The change in 
static stability with latitude is a minimum in July. From July and 
August to September and October, there is a substantial rise in the 
static stability at both latitudes and a sharp increase in the slope of 
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SOUTH LATITUDE 


Fig. 5.5 Variation of static stability slope with latitude in the 
layer 50-100 mb in the Southern Hemisphere from June to 
October 1957. 



Fig. 5.6 Static stability in the layer 50-100 mb as a function of 
latitude in the Southern Hemisphere on 15 October 1957. 
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the profile. In September there is a 15% decrease in static stability 
between 70°S and 90°S. The increased slope in Figure 5.5 tends to 
produce a potential vorticity profile which meets the derived instability 
criterion. Therefore, although in midwinter the flatness of the static 
stability profile in high latitudes prevents the instability criterion 
from being met, in spring, when the vortex is observed to break down, 
a steep slope with latitude develops so that the instability criterion 
tends to be met. This characteristic is also evident in Figure 5.6 
which is produced from a cross-section presented by Palmer and Taylor 
(1960) for conditions in the Antarctic stratosphere just prior to a 
venial breakdown. The decrease with latitude from 70° S to 90° S is again 
obvious and approximately equal in magnitude to the slope given by 
Mahlman (1966) for a Northern Hemisphere, midwinter breakdown. 

5.1.4 Conclusion 

The midwinter stratospheric vortices of the Northern and Southern 
Hemispheres have been shown to have different structures so that the 
Northern Hemisphere vortex tends to meet a previously-derived instability 
criterion, while its southern counterpart does not. This may be the 
reason for the observed difference in breakdown climatologies of the 
two hemispheres; but, a complete explanation must take into account the 
propagation of energy from the troposphere to the stratosphere observed 
during Northern Hemisphere breakdowns (see Matsuno, 1971). The change 
from winter to spring in the static stability distribution in the 
Southern Hemisphere agrees with the observed non-breakdown condition in 
midwinter and the occurrence of the final collapse of the southern 
vortex in spring. 
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5.2 Mean meridional circulation in the wintertime stratosphere 

5.2.1 Background 

As discussed in Section 2.6, very little is known concerning the 
mean meridional circulation of the Southern Hemisphere stratosphere. It 
has been postulated on the basis of surface ozone measurements over 
Antarctica (see Reiter, 1969, 1971) that there is a general sinking over 
the pole during winter. The mean meridional motion pattern in the 
Northern Hemisphere winter has been studied by a number of authors (e. g., 
Vincent, 1968 and Mahlman, 1966). Their studies all reveal a two-celled 
structure with rising motion over the pole and in low latitudes, and 
sinking motion in mid- latitudes. In this section the SIRS-based atmo- 
spheric structure in the stratosphere is used to make computations of 
the zonally-averaged vertical motion fields in both hemispheres during 
a winter month. In the Northern Hemisphere the calculations are also 
carried out using conventional data. 

5.2.2 Method of calculation 

The technique used here to calculate the zonally-averaged vertical 
motion is similar to that used by Mahlman (1966) for the Northern Hemi- 
sphere with conventional data. It is basically a heat budget technique. 
The total derivative of potential temperature (6) is expanded to give 


d0 

at 


30 

3t 


= H + V *V0 + w 


30 

3p 


5.1 


where W is vector horizontal wind on a surface of constant pressure, 
P 

and <i> is the vertical velocity in pressure coordinates. Expansion of 
the advection term in equation 5,1 gives 
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Rearrangement of the terms gives 


5.2 
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Next, equation 5.3 is integrated over the area bounded by two latitudes 
<j>l and 412 , which in this study are usually 10° apart. This results in 
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where a is the Earth’s radius and [ implies averaging around a lati- 
tude circle. The notation used here is explained in detail in Chapter 
4 (Section 4. 1.1) . 

The averaged quantities in the last two terms of equation 5.4 are 
expanded to give 


[V8) x = [v] x [6] x * [Cv) x (6) x ] x 
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Substituting equations 5.5 and 5,6 into equation 5.4, and using 
equation 5.7 results in 



<• i [% 2 ) x (e) x ] x - cos+j 

* A jcos* 2 (v^ 2 H8] x - cos*! [v 4 , 1 l x [e] x J 

5.8 
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In equation 5.8, the second and third lines represent the effect of the 

divergence of the eddy flux and of the mean motion flux, respectively. 

The fourth and fifth lines represent the effect of divergence. 

The last term in equation S.8 involves the correlation of 6 and 

V-V . It is assumed small and, therefore, neglected. The second-to- 
P 

last term in equation 5.8 is well approximated by 
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where [6]. , is the average potential temperature for the area between 

A , (p 

((>1 and $ 2 * The right side of equation 5.9 is integrated to obtain 
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Using the right side of equation 5.10 to replace the second-to- 
last term gives 
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The last two lines in equation 5.11 are of opposite sign and are nearly, 
but not quite, equal. With a temperature gradient of 10°K/10° latitude 
and with 6 approximately equal to 200°K, [0]^ ^ is different from either 
0 at 4>i or 8 at $2 only by about 2 or 3%. The last two lines in equation 
5.11 will thus tend to cancel and their sum is negligible. Therefore, 
they are neglected. 

The integrals in equation 5.11 are evaluated, for example, by the 
approximation 
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where 



is evaluated at latitude 


examp led by equation 5.12, equation 



* ^ 2 . With the approximations 
5.11 becomes 


{cos* 2 [(v +2 ) x ( 8) x ] x - 
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5.13 


The zonally- averaged vertical motion, [w]^, is calculated from equation 
5.13 assuming that w and j^- |~| are uncorrelated. The first term on 
the right side is the local temperature change, and the second term is 
the diabatic heating. The third term is the temperature change due to 
the divergence of the horizontal eddy heat flux. 

The local change in temperature, the static stability and the 
horizontal eddy heat fluxes (geostrophic) can be determined from the 
SIRS-derived atmospheric structure. The diabatic heating, is due 
to radiative heating or cooling and cannot be calculated without a 
knowledge of the vertical distribution of ozone. However, stratospheric, 
radiative heating rates, calculated from climatological conditions, are 
available. Rodgers (1967) has calculated radiative heating rates in 
the troposphere and stratosphere in the Northern Hemisphere up to 70°N, 
for the months of January, April, July and October. Similar calculations 
have been carried out by Dopplick (1972) between 80°N and 80 S for each 
month of the year with the results presented as cross-sections of 
seasonal-average values. The values of Rodgers for January appear to 
be almost identical to those of Dopplick for the Decemb er- January - 
February period in the Northern Hemisphere. Table 5.1 gives the values 
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of || (adapted from Dopplick, 1972) used in this study as a function of 
latitude and layer. The Northern Hemisphere heating rates are those 
given by Dopplick for the December through February period, while the 
June through August values are used for the Southern Hemisphere. 



Latitude 


80°N 

70°N 

60°N 

50°N 

* ’ " 

40°N 

10-30 mb 

-1.2 

-1.2 

-1.0 

-0.8 

-0.5 


-0.9 

-0.9 

-0.8 

-0.6 

-0.3 

50-100 mb 

-0.6 

-0.6 

-0.5 

-0.4 

-0.2 








80°S 

70°S 

60°S 

50°S 

40° S 


-0.6 

-0.8 

, - 

-0.5 

-0.3 

-0.2 

30-50 mb 

-0.4 

-0.4 

-0.3 

-0.2 

-0.2 

SMI 

■M 

-0.3 

-0.3 

-0.3 

1 

-0.2 

-0.2 


Table 5.1 Values of -r— used in vertical motion calculations 
(Adapted from Dopplick, 1972). Units:. °C/day. 

The calculations of [to] ^ are made for each ten degrees of latitude 
between 40° and 80° latitude using fluxes at the intervening latitudes 
of 35°, 45°, etc. The computation of u> at 80° uses the flux at 75° and 
the fact that the flux is zero at 90°. An additional calculation is 
carried out at 85° using the eddy flux at 80°, 

5.2.3 Vertical motion patterns 

The local change term in equation 5.13 is of the order of 0.1°K/day, 
while the heating due to the eddy flux convergence is of the order of 
l°K/day in polar latitudes. The term is negative everywhere during 
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the polar night in the stratosphere and is of the order of l°K/day. 
Because the diabatic heating term varies less with latitude than the 
flux divergence term, the latter is most important in determining the 
latitudinal variability of w. For this reason the enthalpy flux as a 
function of latitude is examined. 

Figure 5.7 shows the poleward heat flux in the 10-30 mb layer 
averaged over January 1970 for the Northern Hemisphere and over July 
1969 for the Southern Hemisphere. Both SIRS and NMC curves are presen- 
ted for the Northern Hemisphere. The geostrophic approximation is used . 
for these computations. In both hemispheres the peak poleward flux is 
reached at 60°. It should be remembered that the polar vortex in the 
Northern Hemisphere broke down in early January that year. Consequently, 
the Northern Hemisphere curves may not be representative of normal (non- 
breakdown) conditions. The two SIRS curves peak at almost the same 
value. The total flux convergence between 60° and the pole is approxi- 
mately equal in the two hemispheres. However, the shape of the two 
curves is different. The Southern Hemisphere curve shows the flux 
decreasing very rapidly just poleward of 60°S and gradually levelling 
out. The greatest flux convergence is in the area between 60°S and 75°S. 
The Northern Hemisphere curve has a relatively constant decrease with 
latitude, indicating that the flux convergence is spread out more 
uniformly than that in the Southern Hemisphere. This difference in 
shape of the enthalpy flux curves is also evident in the 30-50 mb and 
50-100 mb layers, but it is not so obvious. 

The enthalpy flux distributions for all three layers for the 
Southern Hemisphere July 1969 are given in Figure 5.8. In each layer 
there is a rapid decrease in the flux poleward of the peak, then a 
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Fig. 5.7 Poleward heat flux in the 10-30 mb layer-January 1970 in the 
N. H., July 1969 in the S. H. 
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JULY 1969 (ENTIRE MONTH) S. H. 
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Fig. 5.8 Poleward heat fluxes in the stratosphere of the S. H. 
for July 1969. 
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flattening out of the curve in very high latitudes. The value at 80 S 
in the 50-100 mb layer is slightly negative, meaning that between 80°S 
and the pole there is a flux divergence. The flux in the 30-50 mb layer 
at 80°S has a very small positive value so that in this layer the flux 
convergence between 80°S and 90°S is also very small. In the Northern 
Hemisphere there is substantial flux convergence between 80°N and the 
pole (see Figure 5.7). Thus there is a significant difference between 
the hemispheres in the patterns of heat flux and heat flux convergence 
in the polar-night stratosphere. 

The results of evaluating the terms in equation 5.13 and solving 
for u in the Northern Hemisphere for January 1970 with both NMC and 
SIRS-based structure are shown in Figures 5.9a and b. The NMC-based 
calculation (Figure 5.9a) shows rising motion in high latitudes and 
sinking in lower latitudes. The axis of ascent slopes poleward with 
height and the maximum upward motion is located at 20 mb at the pole. 

The pattern of rising motion over the pole and sinking in midlatitudes 
is typical of that found by other authors (e.g., Mahlman, 1966) for the 
Northern Hemisphere. Mahlman (1966) also showed that this distribution 
of ascent and descent occurs before, during and after breakdowns. 

Miller, Brown and Campana (1972) using the quasi-geostrophic, adiabatic 
omega equation obtain a cross-section of w for the period 1-15 January 
1970, which also indicates the same general pattern of ascent and 
descent as in Figure 5.9a. Thus the use of equation 5.13 and NMC 
structure appears to produce a cross-section of vertical motion in the 
stratosphere, which closely resembles previous estimates made for 
different years, and also resembles a calculation made for the same year. 
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Figure 5.9b shows the vertical motions for the Northern Hemisphere 
calculated using the heat fluxes, local temperature changes and static 
stability derived from the SIRS-based structure. There are some simi- 
larities and some differences between Figures 5.9a and b. In Figure 5.9b 
there is rising motion over the pole in the higher layers and sinking 
motion in midlatitudes. There is also a small area of ascent at 60°N 
at 75 mb. The maximum rising motion is found at or near the pole and 
in the top layer (10-30 mb). However, the volume of upward motion is 
much less than that shown in Figure 5.9a. Magnitudes of maximum rising 
and sinking motions are underestimated. 

The vertical velocities in this section are presented in terms of 
oj and in units of mb/day. The conversion from to to w = dz/dt varies 
with density. For a temperature of -60°C the conversions are as follows 
for the various heights: 

at 20 mb, 1.0 mb/day = 3.6mm/sec 

at 40 mb, 1,0 mb/day = 1.8mm/sec 5.14 

at 75 mb, 1.0 mb/day - l.lmm/sec . 

Therefore, the rising motion of 0.6 mb/day at 20 mb in Figure 5.9a is 

approximately equivalent to 2.2mm/sec. 

The main factor in the difference between Figures S.9a and b is the 
underestimation of eddy heat fluxes in the SIRS-based calculation, evi- 
dent in Figure 5.7. Although the SIRS-based calculation undervalues the 
magnitude of the actual fluxes, it reproduces the shape of the NMC curve 
and the locations of flux divergence and convergence. This is also true 
for the two lower layers (30-50 mb and 50-100 mb). The underestimation 
by the SIRS-based fluxes varies with latitude and between layers, the 
minimum underestimation between 40°N and 75°N being 56%. To at least 
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partially correct for this effect a factor of 1.8 (based arbitrarily on 
the 56% figure) is applied to the SIRS-based flux values and the vertical 
motion computation is redone. The results are shown in Figure 5.9c. 

The general rising motion in high latitudes at high levels. and the 
sinking motion in mid-latitudes are still present. The region of rising 
motion has expanded and indicates the slope with height of the NMC-based 
cross-section. Magnitudes of maximum rising and sinking motion are also 
in closer agreement. The adjustment factor 1.8 as based on the minimum 
underestimation is chosen because the author of this study believes that 
the larger underestimations present, especially in the lower layers, are 
due to the low density of satellite radiance data in the Northern Hemi- 
sphere for January 1970 relative to the density of observations over 
the Southern Hemisphere in July 1969, the Southern Hemisphere winter 
month to be examined. 

The results of the vertical motion calculations for the Southern 
Hemisphere for July 1969 are shown in Figures 5.10a and b. Part a of 
the figure is based on fluxes calculated directly from the SIRS-based 
structure, and part b is based on calculations made with the fluxes 
adjusted by the factor of 1.8. Both parts of Figure 5.10 show an area 
of descending motion centered at 50°S, an area of ascent centered at 
about 70°S and a small area of descent over the South Pole, below about 
50 mb. With the assumption of ascent in low latitudes (<30°S) , Figure 
5.10 indicates a possible three-celled meridional circulation in the 
lowest part of the stratosphere, and a two-celled structure in the 
middle stratosphere. The Southern Hemisphere pattern is significantly 
different from that in the other hemisphere, particularly with regard 
to the location and strength of the upward branch. While in the 
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Northern Hemisphere the axis of rising motion tilts poleward in the 
vertical, at southern latitudes the axis is approximately vertical. The 
maximum rate of ascent is over the pole in the Northern Hemisphere, but 
is located at 70°S in the other hemisphere.' The magnitude of the maxi- 
mum also seems to be larger in the Southern Hemisphere, 

The boundary in Figure 5.10 between the main areas of ascent and 
descent (60°S) is the approximate latitude of the polar night jet stream. 
The same is true in the Northern Hemisphere. Because the temperature 
decreases toward the pole, relatively cold air is ascending and rela- 
tively warm air is descending. Therefore, there is a conversion of 
zonal kinetic energy (KZ) to zonal available potential energy (AZ) . It 
seems apparent then, that the stratospheric polar jet is maintained by 
a conversion from eddy kinetic energy (KE) to KZ. Therefore, although 
mean monthly maps show a nearly zonal pattern, the eddies must still be 
very significant in the maintenance of the Southern Hemisphere polar 
night jet stream. 

From Figure 5.10 there does not appear to be a general sinking 
motion over the South Pole, except in the lowest part of the strato- 
sphere. However, since Rubin and Vfeyant (1963) have already noted a 
general subsidence below 75 mb to the surface. over the South Pole, it 
is plausible that ozone-rich air in the lower stratosphere is trans- 
ported downward to the surface to produce the midwinter peaks in 
surface ozone observed over Antarctica (Wisse and Meerburg, 1969). 
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5.3 Variations in the circulation of the Southern Hemisphere 

stratosphere in relation to the presence or absence of a minor 
midwinter warming 

S.3.1 Background 

During July 1969, there were significant variations in the circula- 
tion of the Southern Hemisphere stratosphere. These changes were 
related to the occurrence, early in the 'month, of a minor midwinter 
warming identified by Miller, Finger and Gelman (1970) through the use 
of satellite radiance data. Specifically, Miller et al. used the mean 
SIRS channel 8 radiance around the 60°S latitude circle as a gross 
indicator of variations of stratospheric temperature at that latitude. 
This is reasonable since the weighting function (see Figure 3.1) for 
this channel shows nearly all the contribution to the radiance coming 
from the stratosphere. Large increases in Miller's index are associated 
with a contracting of the polar vortex related to a warming event. The 
plot of this indicator versus time for the Southern Hemisphere winter 
of 1969 shows a sharp rise in mean radiance at the very end of June with 
a relative maximum being reached about 10 July. A very sharp decrease 
follows, after which there is a period of low and relatively constant 
mean radiance. At the very end of July, about the 29th, the mean 
radiance again increases sharply. The sharp rises early and late in 
the month are evidence of minor midwinter warmings. The period of low, 
relatively constant values is what will be referred to here as a quiet 
period. 

In this section the stratospheric circulation is compared between 
the early July minor midwinter warming period and the following quiet 
period. Variations of the mean meridional motion, the eddy heat flux 
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and the wave characteristics are examined. Time- longitude diagrams of 
the satellite-derived thickness and height fields are used to study the 
wave activity. 

5.3.2 Wave characteristics 

Conventional data above 100 mb in the Southern Hemisphere are 
generally inadequate to determine many important characteristics of 
wave motions other than those which can be obtained by single-station 
analysis. An exception is an analysis by Phillpot (1969). A few other 
studies have been done using satellite information. Lovill (1972) and 
Reiter and Lovill (1973) show time-longitude diagrams of total ozone as 
determined from the Nimbus 3 Infrared Interferometer Spectrometer (IRIS) 
experiment. Their diagrams indicate rapid eastward movement of waves 
during the Southern Hemisphere winter. Time-longitude diagrams of 
satellite radiance at 60°S in SIRS channels 7 and 8 are given by Miller, 
Finger and Gelman (1970) covering the same period to be discussed here, 
and Fritz (1970) examines the time- longitude distributions of channel 8 
in May 1969. These diagrams of the radiance emanating from the strato- 
sphere also indicate eastward-moving waves. The radiances represent 
a vertical integration of the temperature structure. Because the waves 
present in the stratosphere probably have a phase shift with height, the 
vertical integration inherent in the satellite radiances may mask some 
waves, especially smaller ones, or may misrepresent the actual amplitudes 
at a particular level. For this reason the thermal structure as obtained 
from the satellite radiances by the regression technique should better 
define the wave activity than merely the radiances themselves. In 
addition, the retrieval technique produces mass distribution information 
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(geopotential height fields) so that phase relations between temperature 
and height can be examined. 

Time- longitude diagrams of the 30-50 mb thickness and the 30 mb 
height averaged between 50°S and TO^S are given in Figures 5,11 and 5.12 
respectively. The analysis on the thickness chart uses an interval of 
50 meters, equivalent in this layer to a temperature interval of about 
3.3°C. Most main features are evident on both diagrams; however, a few 
are clearly discernible only on one diagram. The period covered is 
1-31 July 1969. The latitude belt 50°-70°S covers the zone of the 
Antarctic stratospheric jet stream. 

The period 1-10 July is associated with a minor warming. Figures 
5.11 and 5.12 show the period dominated by a stationary wave number 1 
pattern with a warm ridge centered about 200° E. Superimposed upon the 
stationary pattern is an approximate wave number 3 pattern moving 
rapidly toward the east with a phase speed of about 30° longitude /day. 
Comparing the phases of both the stationary wave number 1 and the 
transient wave number 3 on both diagrams reveals that the temperature 
field lags behind the height field. This results in a westward slope 
with height of the waves. Actually the phase difference between the 
height and temperature (thickness) fields is underestimated. The height 
field in Figure 5.12 is at 30 mb; the thickness field (Figure 5.11) is 
for the layer 30-50 mb, thus representing the temperature distribution 
slightly below 30 mb. The 30 mb temperature field would be shifted 
westward slightly from the 30-50 mb thickness field because of the 
westward phase shift with height. Thus the phase lag between the height 
and temperature fields is somewhat larger than indicated by the two 
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Fig. 5.11 Time- longitude diagram of 30-50 mb thickness averaged from 50°-70 <> S in July 1969, Slopes of 
straight lines indicate phase speeds. 
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diagrams. An eddy heat flux towards the pole is associated with this 
temperature-height phase relation. 

On about 9 July the dominant wave number 1 pattern associated with 
the warming begins to dissolve. For approximately the next IS- 18 days, 
the time- longitude diagrams are dominated by eastward moving waves, 
primarily wave number 2, The approximate phase speed of these waves is 
10° longitude/day. There is also evidence in the thermal pattern 
(Figure S.ll) for a shorter wave (wave number = 3) in the period 17-20 
July. The phase speed of this wave is about 40® longitude/day. During 
this time of the month (11th to the 27th), the radiance index of Miller, 
et al. (1970) is low, which is indicative of a quiet period (i.e., no 
warming event) . 

At almost the end of the month Miller’s radiance parameter again 
indicates a warming at 60®S. Accompanying this increase is an apparent 
change in the wave characteristics with a return to a stationary wave 
number 1 pattern with a wave ridge at about 180®-200°E. There also 
appears to be a transient wave, most evident in Figure 5.12, super- 
imposed on the stationary (k * 1) pattern. 

Therefore, at both the beginning and end of the month, a stationary 
wave number 1 pattern is present and is associated with a midwinter 
minor warming. In the middle of the month, during a quiet period, 
transient wave number 2 dominates. Whether this relation is common to 
the entire winter and/or to other winters is unknown. Phillpot (1969) 
in a time-longitude diagram of 30 mb heights at 65°S for late winter 
and spring of 1967 indicates a nearly stationary wave number 1 pattern 
of varying intensity with the ridge at 170°E to 220®E with transient 
waves (k = 2) superimposed. The location of the wave number 1 ridge is 
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approximately at the same longitude where it is found in Figure 5.12. 
Therefore;, it seems that when wave number 1 dominates the stratospheric 
polar vortex, the highest heights in the midstratosphere are located in 
the Pacific sector. 

The waves shown in Figures 5.11 and 5.12 propagate in a medium with 
large vertical and horizontal shears so that the simple Rossby wave 
equation cannot be expected to hold. However, the waves do exhibit 
Rossby-wave-like characteristics. Principally, the phase speed increases 
with wave number. However, there is difficulty in even attempting to 
apply the Rossby wave formula. A suitable zonal wind speed is difficult 
to choose because of large variations in the horizontal and especially 
in the vertical. The mean zonal geostrophic wind speed at 30 mb for the 
entire month at 60°S is 66 m sec" 1 . At 50 mb the speed is lower, less 
than 50 m sec -1 . Using the 30 mb zonal wind speed, the Rossby wave 
equation gives a k = 1,3 for a stationary wave, reasonably close to the 
k = 1 observed. However, the formula gives considerably higher than 
observed phase speeds foT k = 2. Obviously there are other effects 
present other than those incorporated into the Rossby wave formula. 

Using the Petterssen wave equation (Petterssen, 1956, p. 153), which 
takes into account the horizontal wind shear, results in a reduction 
of the calculated phase speed. However, these speeds are still consid- 
erably higher than the observed. 

5.3.3 Mean meridional motion and eddy heat fluxes 

To inspect the variability of the mean meridional circulation in 
relation to the presence or absence of a midwinter minor warming, calcu- 
lations of zonally-averaged to were carried out for two ten-day periods 
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using the method described in Section 5.2.2. The first interval, 1-10 
July, covers the time of the warming event. The quiet period is 
represented by the time span of 18-27 July. The vertical motion patterns 
during the minor warming and the quiet period exhibit significant 
differences. The vertical component of the mean meridional circulation 
during the minor warming is given in Figures 5.15a and b. The pattern 
for 18-27 July 1969 (quiet period) is given in Figures 5.14a and b. 

During 1-10 July, the area of descent over the South Pole below 50 mb 
disappears, and the ascending and descending branches at about 70°S and 
50°S respectively have larger magnitudes than the monthly averages in 
those locations. The ascending branch is still centered between 70°S 
and 75° S. 

During the quiet period (18-27 July 1969), the vertical motion 
pattern shows descent over the pole at all levels of this calculation 
(Figure 5.14a and b) . According to Figure 5.14b, the descending motion 
covered the area poleward of about 80°S with maximum downward velocity 
occurring in the lowest layer of calculation. Consequently, it appears 
that during this time period the Southern Hemisphere stratosphere has a 
three-celled mean meridional circulation pattern, assuming an area of 
ascent in low latitudes. The ascending branch at 70°S and the descending 
branch farther equatorward are less intense than the monthly average 
(Figure 5.10) and much less intense than they are during the minor 
warming (Figure 5.13), 

The difference in vertical motion patterns between the periods 
1-10 July and 18-27 July 1969 is closely related to the eddy heat flux 
distributions for those two periods (see equation 5.13). These distri- 
butions are shown in Figures 5.15a and b. The peak fluxes during 
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1-10 July are approximately three times as great as those for the 18-27 
July period. Thus the flux divergences on either side of the peak are 
also about three times as great. In very high latitudes another impor- 
tant variation is that, although both the 50-100 mb and 30-50 mb layers 
exhibit eddy flux divergence between SO^S and 90°S during the period 
18-27 July, all layers exhibit flux convergence in that location for 
the 1-10 July period. This is the reason for the existence of the 
difference in vertical motion patterns over the pole between the two 
periods. 

The variations with time of the heat flux and vertical motion 
patterns just discussed lead to some speculation concerning stratospheric 
events leading to the observed mid-winter peaks in surface ozone observed 
over Antarctica (Wisse and Meerburg, 1969). During a minor mid-winter 
warming the poleward heat flux is large and there is apparently no 
descending motion over the South Pole. Assuming that large heat fluxes 
are coincident with large poleward fluxes of ozone in the stratosphere, 
one may postulate a build-up of ozone in the lower stratosphere over 
Antarctica. During an immediately subsequent quiet period the built-up 
ozone reservoir could be tapped by the descending motion over the South 
Pole which has been shown to extend through the troposphere to the 
surface (Rubin and Weyant, 1963) . 

Fritz and Soules (1970) examine a number of Southern Hemisphere 
stratospheric temperature fluctuations during the fall, winter and spring 
of 1969 using SIRS channel 8 data. Between 25 June and 10 July, they 
indicate a large increase in channel 8 radiance in the southern polar 
regions. Accompanying the warming (increasing radiance) in the high 
latitudes of the Southern Hemisphere, they also find decreasing radiance 
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(cooling) in Southern Hemisphere low latitudes and in the Northern 
Hemisphere. The early July case is just one of several similar compen- 
sating changes noted by Fritz and Soules. Another rapid increase of 
zonally-averaged radiance occurs between 26 July and 8 August in high 
southern latitudes. These two periods of radiance increase noted by 
Fritz and Soules have already been discussed in the previous section 
with regard to wave characteristics. In this section it has been noted 
that there are significant differences in the eddy flux of enthalpy 
patterns and in the derived mean meridional motion patterns between the 
time of the early July warming (1-10 July) and the quiet period (18-27 
July) . These observations support the Fritz and Soules postulate that 
irregular (non-seasonal) fluctuations in stratospheric temperature are 
related to changes in large-scale eddies and/or changes in the mean 
meridional circulations. 

5.4 Chapter summary 

The most significant results of this chapter are those regarding 
the structure in both hemispheres of the stratospheric polar vortex and 
the calculation of the mean meridional circulation in the Southern 
Hemisphere. The temperature structure of the southern vortex is shown 
to be significantly different than that in the Northern Hemisphere in 

i 

midwinter so that a previously- derived instability criterion tends not 
to be met in the Southern Hemisphere. This agrees with the absence of 
complete, midwinter breakdowns south of the equator. The Southern 
Hemisphere vortex structure is also shown to vary between midwinter and 
spring so that the instability criterion tends to be met in spring when 
the southern vortex breaks down. . 
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The calculation of zonal ly-averaged vertical motion in the Southern 
Hemisphere stratosphere presented in Section 5.2.3 is the first realistic 
(that is, including the effects of eddy transports) estimate of the mean 
meridional circulation at those heights south of the equator. The 
results indicate the presence of a two-celled structure for the monthly 
mean pattern with rising motion centered about 70°S and the axis of 
sinking motion at approximately 50°S. With temperature decreasing 
toward the pole in these latitudes of the stratosphere the circulation 
is indirect. Consequently, the Antarctic stratospheric jet stream at 
60°S is not maintained by the mean motions, but by the large-scale 
eddies. 

In Section 5.3 there are shown to be significant variations in the 
Southern Hemisphere circulation during July 1969 in relation to the 
presence or absence of a minor midwinter warming. During the warming a 
stationary wave number 1 pattern dominates and the poleward eddy heat 
flux is large, while during a quiet period transient waves of wave number 
2 are most evident and the heat flux is relatively small. The mean 
meridional circulation during the warming (Section 5.3.3) is similar to 
the monthly-mean pattern (Section 5.2.3); but during the quiet period 
the mean circulation becomes, three-celled with the addition cell occur- 
ring in very high latitudes with sinking motion over the South Pole. 

It is possible that the quiet period may be more representative of 
"normal” conditions in the Southern Hemisphere stratosphere than the 
monthly-mean calculation. Therefore, the final determination of the 
normal configuration of the mean meridional circulation will have to 
await additional study. 



CHAPTER 6 


SUMMARY AND CONCLUSIONS 

The basic purpose of this study is to compare the structure and 
flow characteristics of the upper troposphere and stratosphere of the 
Northern and Southern Hemispheres. Atmospheric structure obtained from 
satellite, multi-channel radiance data has provided the primary basis 
for the comparisons. The results of this study confirm some previously- 
made conclusions and assert new conclusions regarding inter-hemispheric 
comparisons and the Southern Hemisphere circulation. 

The satellite, multi-channel radiance data have been shown to be 
a valuable asset in the study of the general circulation. In this 
study they have provided information on atmospheric structure with 
equal, or nearly equal, time and space resolutions in the two hemi- 
spheres. This investigation is, therefore, the first comparison of 
circulations of the Northern and Southern Hemispheres based on a data 
system having comparable resolution and quality of data above and below 
the equator. The radiosonde network of the Southern Hemisphere has 
nearly an order of magnitude fewer stations than the Northern Hemisphere 
version. Therefore, the comparison between the hemispheres of statis- 
tics based on these radiosonde systems may be open to question. This 
may be especially true of quantities whose estimation requires good 
horizontal resolution on the time scale of a day. 
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Most previous observational studies of the atmosphere’s general 

circulation have centered on the Northern Hemisphere, because of the 

) 

better upper-air network there. However, 'in order to understand the 
general circulation on a global basis we must first develop a knowledge 
of the characteristics of the circulation in both hemispheres. Because 
there is little hope that the number of Southern Hemisphere radiosonde 
stations will increase to become comparable to the network in the 
Northern Hemisphere, there will be continued heavy reliance on satellite 
data to study the general circulation below the equator. Numerical 
models of the global general circulation in the troposphere and strato- 
sphere now exist. In the Northern Hemisphere these models can be 
verified against statistics compiled from the relatively dense radio- 
sonde network; in the Southern Hemisphere this verification will have to 
depend more heavily on satellite-based studies, especially for eddy 
terms, and especially for layers above the mid-troposphere . Hopefully, 
the current investigation has provided a beginning in that direction. 

The results of this investigation do show that the satellite 
radiance data, as used here, have two main limitations in this type of 
study. First, because of cloud effects in the lower levels, results 
for the most part are restricted to layers above 500 mb. Even above 
500 mb results improve with height into the stratosphere. The second 
limitation is the tendency of the satellite-derived atmospheric structure 
to underestimate the amplitude of eddies. 

The major results and conclusions of this study are summarized in 
the following paragraphs. 

Zonal ly-averaged vertical motion fields in the stratosphere calcu- 
lated in both hemispheres confirm a two-cell structure to the mean 
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meridional circulation in Northern Hemisphere winter with rising motion 
over the North Pole. In the Southern Hemisphere winter, the calcula- 
tions also show a two-cell structure in the mid- stratosphere, but with 
the main ascending branch located at approximately 70° S. The mean 
circulation in the vicinity of the Antarctic stratospheric jet stream 
is indirect, thus the maintenance of the jet stream must depend on the 
large-scale eddies. The calculation of the mean meridional circulation 
of the Southern Hemisphere stratosphere presented in this paper is the 
first realistic estimation of that circulation. That is, it is the 
first to include the very important influence of large-scale eddies. 

There are significant variations during the Southern Hemisphere 
winter in relation to stratospheric synoptic events. During a minor 

t 

midwinter warming, the flow is dominated by a stationary wave number 1 
pattern, the eddy flux of heat at 60 °S is large, and the vertical 
motion magnitudes are also large. During a "quiet" period, the flow is 
dominated by a transient wave number 2 pattern, the eddy flux is 
relatively small, and a descending branch of the mean meridional circu- 
lation appears over the South Pole. Consequently, during this quiet 
period there appears to be three instead of two cells in the mean 
meridional circulation. 

A previous calculation based on IGY data in the Southern Hemisphere 
concluded that the Southern Hemisphere mid- latitude eddy flux of heat 
is approximately one half as large as that in the Northern Hemisphere. 
Another study, based on energy budget considerations, concludes that 
the two values are about the same. The results of the present study 
indicate that the eddy flux is much closer to being equal in the two 
hemispheres than to being different by a factor of two. 
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The existence of longitudinal variations in the pattern of eddy 
heat flux in southern mid- latitudes , as there are in northern mid- 
latitudes, is also shown. These variations are best defined in both 
hemispheres during winter. The Northern Hemisphere variations are 
closely tied to standing waves; in the Southern Hemisphere the varia- 
tions are due more to longitudinal variations of the strength of 
transient eddies. In winter, at 45°S, three peaks are well defined, 
seemingly associated with the three land areas, Australia, South 
America and Africa. 

The relative importance of standing eddies in the two hemispheres 
has been examined closely in this study. The calculations confirm that, 
on a hemispheric-averaged basis, the Southern Hemisphere has a lower 
percentage of AE and KE contained in standing eddies. In the Southern 
Hemisphere the percentage of energy contained in the standing eddies is 
smaller in summer than in winter; in the Northern Hemisphere, the 
percentage either remains the same or increases from winter to summer, 
depending on the layer and parameter under consideration. This differ- 
ence reflects the dominance of the Asian summer monsoon in the Northern 
Hemisphere. 

The distribution with latitude in the upper troposphere of the 
ratio of standing eddy temperature variance to the total variance shows 
large inter-hemispheric differences. This parameter indicates the 
importance of standing eddies as a function of latitude. In winter, 
poleward of 30°, the Northern Hemisphere shows one peak at about 5Q°N; 
the Southern Hemisphere has two peaks, one at 40°S and one over 
Antarctica. In mid-latitudes (40-60°), the Southern Hemisphere value is 
lower than that in the Northern Hemisphere; in high latitudes (>60°), 
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the reverse is true. In summer, in low latitudes (<35°), both hemi- 
spheres show maxima associated with the summer monsoon in each 
hemisphere. The Northern Hemisphere maximum is much larger than the 
Southern Hemisphere one. At all latitudes in summer the ratio is 
larger in the Northern Hemisphere. 

In winter, at the latitude of maximum eddy kinetic energy (~30°), 
the standing eddies contribute substantially to the total. In the 
Southern Hemisphere, the percentage contribution is about 30%; in the 
Northern Hemisphere it is about 40%. These percentages are larger than 
the hemispheric-averaged values. In summer the percentages are less 
in both hemispheres . 

In the Northern Hemisphere during winter, about 50% of the eddy 
heat transport in mid-latitudes is accomplished by the standing eddies. 
In the Southern Hemisphere, the percentage contribution is much less. 
According to calculations in this study, the Southern Hemisphere per- 
centage is about 20%. Other previous studies have indicated a lower 
percentage. 

The satellite-based zonal mean temperature and mass structure 
confirm the inter-hemispheric differences previously noted from 
climatological information. The upper tropospheric temperature gradient 
between 20° and 80° latitude is larger in the Southern Hemisphere in 
both winter and summer, with the largest disparity occurring in summer. 
The same difference is evident in the distribution of 200 mb geopoten- 
tial height. These variations are also reflected in the computed values 
of zonal available potential energy (AZ) and zonal kinetic energy (KZ) , 
In summer the AZ of the Southern Hemisphere is twice as large as that 
in the Northern Hemisphere. Because the winter difference is relatively 
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small, there is an annual cycle of global AZ with peak values occurring 
during the northern winter. 

The contributions from the upper tropospheric and tropopause layers 
to the eddy available potential energy (AE) and eddy kinetic energy (KE) 
do not show the large inter-hemispheric differences that are evident in 
the zonal components, AZ and KZ. Therefore, the ratios of AZ to AE and 
KZ to KE are larger in the Southern Hemisphere. 

The winter stratospheres of the Northern and Southern Hemispheres 
are also studied in relation to the difference in breakdown clima- 
tologies. The temperature structure of the winter stratospheric vortices 
are shown to vary between the hemispheres in a way that a previously- 
derived instability criterion tends to be met in the Northern Hemisphere, 
but not in the Southern Hemisphere, during midwinter. This difference 
may be related to the dissimilar midwinter breakdown climatologies. In 
the Southern Hemisphere, the structure varies between midwinter and 
spring so that the vortex tends to meet the instability criterion in 
spring when the vortex normally breaks down. 
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APPENDIX 


REGRESSION COEFFICIENTS 

The general form of the regression equations used in this study is 

D - ct q + ctjNi + Q 2 N 2 + + ot gN 8 > A. 1 

where D is the geopotential thickness, <*0 is the equation constant, 04 
through otg are the regression coefficients and Nj through Ng are the 
radiances in the eight SIRS channels. 

In the four tables which follow the coefficients in the thickness - 
radiance regressions equations are given. Coefficients are presented 
for each of the eight layers and for each of three latitude bands : 

I, IQ* ) II, 40°-60°; III, 60°-80°. The variance explained (V. E.) 
and the standard error of estimate (S. E. E.) in both meters and "C are 
also exhibited. No regression equations were determined for the 
10-30 mb layer in the Southern Hemisphere, because of a lack of conven- 
tional data at that altitude. 

Each equation's constant, ao, is given in the column marked 
"Constant." The values of 04 through ag are shown in the following 
eight columns. A blank location means that the a for that channel is 
zero. The units in the regression equations are: thickness, m; 

temperature, °C; radiance, ergs cm -2 sec -1 steradian -1 (cm -1 ) . A 

typical value for the radiance N is 50 ergs cm -2 sec -1 steradian -1 
(cm - 1 ) 
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An examination of the coefficients in the tables indicates that 
layers in the stratosphere have regression equations with large positive 
coefficients for channels six and seven. This is physically sound, 
since the weighting functions for those channels have peaks in the lower 
stratosphere (see Figure 3.1). As one moves from stratospheric layers 
to those below 100 mb, the major positive coefficient shifts toward 
channels with weighting functions which peak at lower altitudes. 

However, in the two lowest layers, especially 700-1000 mb, this physical 
relationship does not seem to hold. In other words, based on the dis- 
tribution of weighting functions in Figure 3.1, one would expect a large 
positive coefficient for channel one, two or three in the last row of 
each of the tables. This, however, is not the case. Because not all 
cloud-contaminated data have been eliminated, the selection procedure 
picks other channels which are not as affected by clouds as are channels 
one through three. Since there is a significant negative 
correlation between tropospheric and stratospheric temperature, channel 
five, six or seven enters the regression equation with the dominant 
coefficient, a negative one. 
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Table A.l Regression coefficients for July 1969, Northern Hemisphere, 



— - 

Constant 

— 

— 

n 2 

«3_ 



“6 

— *2 — ■ 

*8 

v,e.c%) 

S.E.E. 

W 

S.E.E. 

fC) 

10- 

I 

6449.2 








13.29 

60.0 

24.8 

0.3 

30 mb 

II 

5918.7 





-54.29 


86 , 86 

-9.49 

76.4 

108.1 

3.4 


m 

7084,6 





-53.69 


S9.06 

-5.81 

66. S 

97.1 

3.0 

30- 

i 

2914.4 







25.61 

-15.33 

64.4 

42,0 

2.8 

50 mb 

n 

2655.6 







29,7$ 

-14.91 

6S.3 

54.0 

3,6 


m 

2540,0 







32.69 

-IS. 69 

PS. 4 

56.9 

2.5 

SO- 

i 

3224.9 





-14,61 


70.12 

-27.81 

94.5 

37.9 

1.9 

10Q mb 

ir 

35 19. 2 






37.25 


-14.79 

76.3 

52.6 

2.6 


ni 

3146,7 






30,62 


13.20 

-13.54 

87.7 

53,9 

2.6 

100- 

... .. . 
i 

3253.7 




-7.27 


47.51 

- 

-10.36 

81.6 

38.7 

1.9 

200 mb 

n 

3554.5 




-16.16 

55,28 


-13,96 

-2,67 

82.5 

40.2 

2.0 


m 

3463.6 




-16.30 

58.57 


-18.64 


89.0 

43,7 

2.2 

200- 


1896.9 


-1.15 


8.00 

2D. 76 


-12.15 


70.3 

1 

26.9 , 

2.3 

300 mb 

11 

1969.0 




-4.00 

28.38 


-11,83 

1.84 

70.7 

27. S 

2.3 


m ! 

2162.2 



-2,80 


21.75 

- . 

-8.55 


72.2 

32.7 , 

2.7 

300- 

i 1 

3677.8 


-1.12 


| 22.06 


-26.53 


i 

79. S 

40.7 

2,7 

500 mb 

ii 1 

2934.3 


-4.22 


| 41.45 


-27.43 



64.0 

45.3 

3.0 


in 

3004.3 


-2 . 56 


i JS.48 

-22.04 


-3.63 


51.6 

38. S 

2,6 

SQO- 

i 1 

2603.5 




10.11 


-27.62 


7.7B 

BO.O 

28.9 

2,9 

700 rab 

n 

2117.9 


-3.32 


35.51 


-34.75 

8. 32 


j 62,9 

40.8 

4,2 


m ' 

i 

2110.0 




21.88 


-26.04 

9,04 


1 SB .2 

38.6 

3.9 

700- 

i i 

3524.7 


2.41 




-26.57 


[ 5.59 

64,1 

44.1 

4,2 

1000 mb 

ii i 

2456.1 

1 



23.90 


-14.48 

-9,29 

3.54 

49, a 

66.0 

6. 3 

k- —— — 

in 

2564.4 

[ 3.36 

— . 

4.95 



-9.10 



L ■-» 

! 58.2 

59.5 

5.7 


Table A. 2 Regression coefficients for January 1970, Northern 
Hemisphere . 
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Table A. 3 Regression coefficients for July, 1969, Southern 
Hemisphere. 



Table A. 4 Regression coefficients for January, 1970, Southern 
Hemisphere. 
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